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Abstract 



In quantum field theory in curved spacetimes the construction of the algebra of observables 
of linear fields is today well understood. However, it remains a non-trivial task to construct 
physically meaningful states on the algebra. For instance, we are in the unsatisfactory 
situation that there exist no examples of states suited to describe local thermal equilibrium 
in a non-stationary spacetime. In this thesis, we construct a class of states for the Klein- 
Gordon field in Robert son- Walker spacetimes, which seem to provide the first example 
of thermal states in a spacetime without time translation symmetry. More precisely, in 
the setting of real, linear, scalar fields in Robertson- Walker spacetimes we define on the 
set of homogeneous, isotropic, quasi-free states a free energy functional that is based on 
the averaged energy density measured by an isotropic observer along his worldline. This 
functional is well defined and lower bounded by a suitable quantum energy inequality. 
Subsequently, we minimize this functional and obtain states that we interpret as 'almost 
equilibrium states'. It turns out that the states of low energy, which were recently introduced 
in [Olb07bj . are the ground states of the almost equilibrium states. Finally, we prove that 
the almost equilibrium states satisfy the Hadamard condition, which qualifies them as 
physically meaningful states. 

Zusammenfassung 

In der Quantenfeldtheorie in gekriimmter Raumzeit ist die Konstruktion der Algebra der 
Observablen linearer Felder heutzutage gut verstanden. Es ist jedoch eine nicht-triviale 
Aufgabe physikalische Zustande auf der Algebra zu konstruieren. Zum Beispiel sind wir 
in der unzufriedenstellenden Situation, dass keine Zustande bekannt sind, die es erlauben 
ein lokales thermales Gleichgewicht in einer nicht-stationaren Raumzeit zu beschreiben. 
In dieser Arbeit konstruieren wir ein Klasse von Zustanden fiir das Klein-Gordon Feld in 
Robertson- Walker Raumzeiten, die das erste Beispiel thermaler Zustande in einer Raumzeit 
ohne Translationssymmetrie bezuglich der Zeit darstellen diirften. Genauer gesagt definie- 
ren wir fiir das reelle lineare skalare Feld in Robert son- Walker Raumzeiten ein Funktional 
fiir die freie Energie auf der Menge der homogenen isotropen quasifreien Zustande, welches 
auf der gemittelten Energiedichte basiert, die ein isotroper Beobachter auf seiner Weltlinie 
misst. Dieses Funktional ist wohldefiniert und von unten beschrankt, dank einer geeigne- 
ten Quanten-Energie-Ungleichung. Danach minimieren wir dieses Funktional und erhalten 
so Zustande, welche wir als Fast-Gleichgewichts-Zustande ("almost equilibrium states") 
interpretieren. Es stellt sich heraus, dass die Zustande niedriger Energie ("states of low 
energy"), welche kurzlich in [Olb07b| definiert wurden, die naturlichen Grundzustande der 
Fast-Gleichgewichts-Zustande sind. SchlieBlich beweisen wir, dass die Fast-Gleichgewichts- 
Zustande die Hadamard-Bedingung erfiillen, was sie als physikalisch sinnvolle Zustande 
auszeichnet. 



Dedicated to the memory of Bernd Kuckert. 
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1 Introduction 



There are plenty of reasons to investigate quantum fields in curved spacetimes. Accepting 
for the time being that we don't know how to formulate a quantum theory of gravity, 
we may still have the desire to investigate the physics in the intermediate regime, where 
the quantum nature of fields is already apparent and spacetime is curved, but the quantum 
nature of gravitation plays no crucial role. That such a theory can produce interesting results 
is undoubted since the discovery of the Hawking effect, which establishes a remarkable 
connection between two seemingly distinct branches of physics, namely, gravitational (black 
hole) physics and thermodynamics. The latter field is also of major importance in modern 
cosmology, where the cosmic microwave background (CMB), its isotropy on large scales, 
and its anisotropy at small scales remains a startling puzzle. It is also remarkable that 
cosmological observations, e.g., the ones leading to the postulation of dark matter and dark 
energy, nowadays become a major trigger for particle physics research. For example, the 
Large Hadron Collider (LHC) at CERN, scheduled to start operating this year, is hoped to 
provide clues on the constituents of dark matter. Surely, the interplay between cosmology 
and particle physics will increase in the years to come. 

The formulation of quantum field theory in curved spacetime is hampered by the fact 
that the traditional formulation relies heavily on Poincare symmetry, which is lacking in a 
general curved spacetime. Consequently, notions like 'particles' and 'vacuum', which de- 
pend on Poincare symmetry, are not well defined in this setting. Furthermore, quantum 
field theory becomes ambiguous because of the appearance of unitarily inequivalent Hilbert 
space representations of the algebra of observables. An expanding, homogeneous, isotropic 
universe is described by a Robertson- Walker metric, which is non-stationary and comes in 
three types, of which two have curved spacelike sections. Owing to the present symmetry, 
the curvature of the spacelike sections does not pose a serious problem for the formulation 
of quantum field theory. On has analogues of the Fourier transform, and the quantum field 
can be written, as usual, as an integral over modes, albeit the coefficients in the mode 
decomposition, which in stationary spacetimes are interpreted as creation and annihilation 
operators, loose their meaning. However, non-stationarity, i.e., the absence of time transla- 
tion symmetry, brings about, among others, the problems mentioned above. In particular, 
the lack of a timelike Killing field results in a failure of energy conservation and, conse- 
quently, the lack of a Hamiltonian, which generates time-evolution in the non- stationary 
case. 

A thorough analysis of these problems is viable within the algebraic approach to quantum 
field theory, which can be rigorously formulated for quantum fields on curved backgrounds. 
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Chapter 1. Introduction 



In this approach, one first constructs the algebra of observables for the quantum system, 
satisfying, for instance, the canonical commutation relations. In a second step, one con- 
structs states with desired properties as linear functionals on the algebra. This course of 
action disentangles, effectively, problems tied to the algebra of observables from problems 
arising at the level of representations and states. While the first step is, at least for free 
fields, to a large degree understood, the second one is, put mildly, less explored. One has 
to admit, that the number of known examples for physically meaningful states in curved 
spacetimes is small. For example, to date, there are no explicitly known thermal states with 
respect to global time evolution in a non- stationary spacetime. 

A different problem - which is already present in flat spacetime -, namely, the lower 
unboundedness of point-like energy densities in quantum field theory has led to the so- 
called quantum inequalities. These uncertainty-type inequalities give lower bounds on the 
averaged energy densities measurable in (a class of) quantum states, where the averaging 
procedure involves a sampling function with suitable properties. Quantum inequalities were 
established in flat spacetimes as well as in arbitrary globally hyperbolic spacetimes. They are 
important in general relativity since they reestablish the macroscopic energy conditions, e.g., 
the weak energy condition, utilized for proving singularity theorems. Furthermore, it has 
been argued that they may serve as stability conditions for quantum systems on a mesoscopic 
scale between the Hadamard condition (microscopic scale) and passivity (macroscopic scale) 
[FV03] . 

The initial motivation for the project presented in this thesis can be stated as follows. 
Since there are, by quantum inequalities, locally meaningful energy-quantities even in curved 
spacetimes, how do the ground states and (thermal) equilibrium states of these 'energies' 
look like? Regarding ground states, the analysis carried out by Olbermann in [ Olb07b] 
in the setting of real, linear, scalar fields in Robertson- Walker spacetimes resulted in the 
definition of states of low energy. These are pure, homogeneous, isotropic, quasi-free states 
that minimize the averaged energy density measured by an observer along a timelike curve. 
Moreover, it could be proven, by comparison with adiabatic vacuum states, that the two- 
point distributions of the states of low energy satisfy the Hadamard condition, which is a 
condition on the short distance behaviour of the two-point distribution to be satisfied by 
physically reasonable states. Encouraged by this result, we construct, in the same setting, 
a class of (thermal) almost equilibrium states with respect to the timelike averaged energy 
density. 

The basic idea of our construction is the following. For a finite quantum statistical system 
in contact with a much larger reservoir at a constant positive temperature 4, the free energy 
T = £ — 4<S, where £ is the inner energy and S denotes the entropy, is the maximum 
amount of work that the system can perform. It turns out that the equilibrium state of 
such system is characterized uniquely by the property that it minimizes the free energy. 
Following this general principle, we define on the set of homogeneous, isotropic, quasi-free 
states of the Klein-Gordon field in Robert son- Walker spacetimes a free energy functional, 
where a timelike averaged energy density takes the role of £. We show that, owing to the 
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symmetry of the spacetimes, this energy quantity has the form of an integral over the modes 
of the quantum field. In view of this, we associate to each mode an entropy functional and 
require the mode to evolve according to the KMS condition, i.e., as a system in thermal 
equilibrium at inverse temperature (3. By a certain quantum inequality proved in |Few00j . 
our free energy functional is bounded below within the class of Hadamard states. Indeed, 
we succeed in calculating the two-point distribution of the unique state that minimizes our 
free energy functional for a given sampling function. We call the resulting state an almost 
equilibrium state a inverse temperature f3 associated to the sampling function used for 
averaging the energy density. We finally prove that the almost equilibrium states satisfy the 
Hadamard condition. The last step is crucial since, as noted above, the quantum inequality 
on which we base our definition is valid only within the class of Hadamard states. At the 
same time, the Hadamard property proves the almost equilibrium states to be physically 
meaningful states. 

Outline The outline of this thesis is as follows. In section 2 we give an account of the basics 
of quantum field theory in curved spacetimes. We consider both the quantization scheme 
via construction of a Hilbert space as well as the algebraic method. We introduce the notion 
of quasi-free states and Hadamard states. The main sources for this chapter are [Wal94j 
and [KW91] but the reader, who is not well acquainted with algebraic quantum field theory 
will also benefit from |Haa96] and [Emc72] . In section 3, we specialize to Robertson- Walker 
spacetimes and recall relevant facts about homogeneous, isotropic, quasi-free states of the 
Klein-Gordon field in this setting. In particular, we introduce the two-point distribution 
of such states. Most of the the material presented here can be found in the seminal paper 
|LR90j . We add to the known facts a characterization of the two-point distribution in 
terms of the mode solutions of the Klein-Gordon equation. Chapter 4 is devoted to quantum 
inequalities. A brief survey of the topic is followed by a quotation of the quantum inequality 
on which our construction is based [FewOOj . Then, we calculate the specific expression for 
the averaged stress-energy tensor, since this will be part of the free energy that we will 
minimize. Section 5 is the main part of this thesis, in which we explain and accomplish the 
construction of the almost equilibrium states. As a necessary step for our construction, but 
maybe also interesting in its own right, we calculate the generator of KMS states on the Weyl 
algebra of a system with one degree of freedom. Then, we define the almost equilibrium 
states and prove that they are indeed Hadamard states. In section 6, a summary of the 
construction is given and supplemented by a few ideas on possible future projects. 
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2 Quantum field theory in curved 
spacetime 



The overall setting of this thesis is quantum field theory in curved spacetime. This is a 
generalization of ordinary quantum field theory in the sense that the fields propagate on a 
curved Lorentzian manifold rather than flat Minkowski space. It may be seen as a semi- 
classical approach to 'quantum gravity' since the fields are quantized but the background 
spacetime is taken to be classical and fixed. Any backreaction of the field on the spacetime 
is, in principle, encoded in the semiclassical Einstein equations 



though, to date, few is known about backreaction effects, due to some serious difficulties in 
calculating (T Mi ,). See the discussion in |Wal94j on this point. 

The absence of symmetries is a big obstacle for the traditional formulation of quantum 
field theory on a general spacetime. However, the algebraic quantum field theory, as defined 
in [HK64j . has a well defined generalization to curved spacetimes |Dim80j . The fundamen- 
tal object in the algebraic formulation is a net of local algebras on a spacetime. Such a 
net associates to each open region O in spacetime an algebra 21(0) of observables that are 
measurable in the region. The latest enhancement of the algebraic approach is the so-called 
locally covariant quantum field theory |BFV03] . In this approach a locally covariant quan- 
tum field theory is defined as a covariant functor between the category of globally hyperbolic 
Lorentzian spacetimes with admissible (isometric, orientation and time-orientation preserv- 
ing, causally convex) embeddings as morphisms and the category of unital C*-algebras with 
certain homomorphisms as morphisms. 

In the first section of this chapter, we introduce the necessary facts about globally hy- 
perbolic spacetimes and the Cauchy problem thereon. Then, we present the quantization 
scheme in curved spacetime, where we follow the presentation found in the monograph 
|Wal94] and the review [KW91j . At last, we introduce the notions of quasi-free states and 
Hadamard states. Our presentation of the general theory is strongly biased by our needs 
in the later chapters. 
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Chapter 2. Quantum field theory in curved spacetime 



2.1 Classical preliminaries 
2.1.1 Spacetime structure 

We adopt the common viewpoint that spacetime is modelled by a smooth (in the sense 
that it is C°°, Hausdorff, paracompact, and connected) four- dimensional manifold Ai with 
Lorentzian metric g^ u . We also assume that the spacetime under consideration is orientable 
and time-orientable. Actually, we will restrict ourselves to the class of Robertson- Walker 
spacetimes, which fulfill these hypotheses. Following the loosely established tradition of 
field theorists, we choose the signature of the metric to be (+, — , — , — ). According to 
this convention, a vector is called timelike if its norm Q flu v fJ- v u is positive, spacelike if the 
norm is negative, and null if it is zero. A curve j(t) : R — > M. with tangent vectors 7 M (t) 
that are timelike everywhere is called a timelike curve. It is called a causal curve if the 
tangent vectors are timelike or null everywhere. For K C M. the causal future J + {K,) and 
the causal past J~{K) are the sets of all points that can be reached from K, by a future/past 
directed causal curve. A point x G M. is said to be the future endpoint of a causal curve if 
for every neighborhood O of x there is a to such that j(t) G O for every t > t . A causal 
curve is said to be future inextendible if it has no future endpoint. Analogously, we can 
define past inextendible curves. A causal curve is called inextendible if it is future and past 
inextendible [HE73] . Less precisely, one can say that an inextendible curve can end only at 
infinity or at some initial or final singularity. The future (+) /past (-) domain of dependence 
-D ± (S) of E is the set of all x G M. such that every past/future inextendible causal curve 
through x intersects E [Wal84j. 

A Cauchy surface for a Lorentzian manifold M. is a spacelike hypersurface which is 
intersected exactly once by every inextendible causal curve in Ai [Ful 89] . Equivalently, 
the domain of dependence -D(E) = Z} + (E) U Z) _ (E) of a Cauchy surface is the entire 
spacetime. A spacetime {M.-,Q^y) is said to be stationary if there exists a one-parameter 
group of isometries E t : M. — > Ai, E^g^ = Q^u, t G R whose orbits are timelike curves. 
Equivalently, there exists a global timelike Killing vector field £ M satisfying Killing's equation 
V M £„ + = 0, namely, the generator of E t . The spacetime is said to be static if it is 

stationary and if, in addition, there exists a spacelike hypersurface E that is orthogonal to 
the orbits of the isometry. In a local coordinate system x\ % = 1,2, 3, the metric of a static 
spacetime may be written as ds 2 = N(x) 2 dt 2 — i)ij(x)dx l dx : ' , where f)jj is the induced metric 
on the Cauchy surfaces, and N G C°°(E) is the lapse function. In this case, the timelike 
Killing vector field satisfies ^V^p] = 0. For completeness, we note that a spacetime is 
called ultrastatic if it possesses a timelike Killing vector field £ M that is orthogonal to the 
spacelike hypersurfaces and obeys j m,£ M £ 1 ' = 1. In this case, the lapse function is N — 1. 

In order to guarantee the well posedness of the classical time evolution of a field, one 
usually assumes globally hyperbolicity of the underlying spacetime. There are different 
equivalent definitions of global hyperbolicity. For example, a spacetime {M.,Q^ U ) is called 
globally hyperbolic if there are no closed causal curves in Ai and the collection of all causal 
curves joining two arbitrary points xi,X2 G A4 is compact (in a suitable topology). Due to 
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a theorem by Geroch |Ger70] an n- dimensional globally hyperbolic spacetime (Ai, Q^) can 
be foliated with a family of Cauchy surfaces that are diffeomorphic to an (n— l)-dimensional 
manifold S. This means that a globally hyperbolic spacetime has the topology Ai = R x E. 
Global hyperbolicity is a condition on the geometry of spacetime that ensures the existence 
and uniqueness of global solutions to hyperbolic differential equations. In a broader sense, 
global hyperbolicity of Lorentzian spacetimes serves as an analog for the notion of com- 
pleteness of Riemannian manifolds [BGP07] . Since quantum dynamics is usually modelled 
around classical dynamics, global hyperbolicity of the underlying spacetime is a reasonable 
assumption in a semiclassical approach like quantum field theory in curved spacetime. We 
refer the interested reader to [HE731 IWal84| IBGP07] for further elaborations on this topic. 

2.1.2 The Cauchy problem 

A thorough and general treatment of the Cauchy problem on globally hyperbolic spacetimes 
can be found in the monograph |BGP07j . The results there are valid for arbitrary (complex) 
vector bundles on general globally hyperbolic manifolds. However, in order to keep in touch 
with the physical literature and to avoid unnecessary complexity, we will use the traditional 
formulation found in [Dim80t IWal94] , keeping in mind the more general results. 

In this thesis, we are concerned with the special case of a real, linear, scalar field on a 
globally hyperbolic spacetime (Ai, Q^u) satisfying the Klein-Gordon equation, 

(□ + m 2 )0 = O, (2.2) 

where D g = g M ^V M V v is the wave operator associated to and m > is the mass 
parameter. By the global hyperbolicity of Ai, the Klein-Gordon equation (I2.2p has a 
well-posed initial value formulation. This means that there exist unique continuous linear 
operators 

E ± ■ 9>{M) -> S(M) , (2.3) 

called the advanced fundamental solution, E + , and retarded fundamental solution, E~ , 
with the following properties: 

(□ B + m 2 )^/ = f = E ± (D s + m 2 )f , (2.4) 
supp(£+/)c J + (supp/), (2.5) 
supp(£T/)c.7-(supp/) (2.6) 

for / E <&(M). We use the notation 9{M) := G^(M) and S\M) := C°°(M), denoting 
the set of (complex-valued) smooth functions (with compact support in the case of Sf) 
on the manifold Ai. The first of these properties says that E + and E~ are solutions 
to the inhomogeneous Klein-Gordon equation. The other properties are sensible support 
properties. Combining the advanced and retarded fundamental solutions, one constructs 
the fundamental solution 



E:=E + -E- : @{M) -> S{M) 



(2.7) 
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of the Cauchy problem, which has the properties 

(□„ + m 2 )Ef = = E(D g + m 2 )f , (2.8) 
supp(£/) C (J+(supp/) U J-(supp/)) (2.9) 

for / G @{M). 

There are different schemes available for the quantization of a classical fields, some of 
which start from complex solutions and some of which use real-valued functions. In this 
chapter we will be concerned with the latter spaces. Denote by 5? C $\M., R) the space 
of classical, real-valued, smooth solutions to the Klein-Gordon equation that have compact 
support on Cauchy surfaces, and denote by £F := R) the test function space of real- 

valued, smooth functions of the fundamental solution, with compact support. Then the 
following lemma states some important properties 

Lemma 2.1 (|Wal94|) The map E : 2? — > 5? satisfies the following three properties 

(i) E is onto, i.e., every G 5? can be expressed as = Ef for some f G . 

(ii) Ef = if and only if f = (D + m 2 )g for some g G . 
(Hi) For all G 5? and all f G , we have 

<f>(f) := J d A x 4>(x)f(x) = a(Ef,4>) . (2.10) 

The fundamental solution plays an important role in the quantization procedures to be 
described later. For example, it determines the values of the commutators of the quantum 
fields. Even more important for us, it can be used to calculate the four-smeared two-point 
distribution from the three-smeared two-point distribution of the quantum fields, as we will 
see in section 13.3.11 

The fundamental solution E maps test functions on a globally hyperbolic spacetime 
to solutions that arise from initial data on some Cauchy surface. More precisely, this 
correspondence can be formulated as follows [Dim80j. Define for a given Cauchy surface £ 
the restriction operator p and the forward normal derivative p\ by 

p : £{M,R) -> <f(E,R) pi : ${M, R) — > <f (£, R) (2.11) 

0^0|e, 0^KV^)| s , (2.12) 

where n M denotes the unit forward normal to £ and n^V^ is the Lie derivative in this 
direction. Now, the following holds. Let £ be any Cauchy surface and let u,p G ^(£, R) be 
any pair of smooth functions with compact support on the Cauchy surface £. Then, there 
exists a unique solution G 5? defined on all of M. to the Klein Gordon equation (12. 2p that 
is related to its initial values on £ by po((f>) = u and pi(0) = p. Furthermore, 

supp C (U± J^suppu)) U (U± ^(suppp)) , (2.13) 

i.e., for any closed subset of £ the solution restricted to the corresponding domain of 
dependence depends only on the initial data in that subset. 
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2.1.3 Phase spaces 

The quantum theory of linear systems is usually modelled after the classical theory as it 
becomes apparent in the 'Poisson bracket goes to commutator' rule. A sophisticated version 
of this rule is used for linear fields in curved spacetimes. Now, we introduce the classical 
part of this correspondence, namely, the classical phase space. 

A classical phase space is a symplectic vector space, i.e., a pair (Y,a), where Y is a 
vector space and cr is a bilinear form a : Y x Y — > R that is symplectic, a(f, g) = —<r(g, f), 
and non-degenerate, which means if a(f, g) = for all / G Y then g = 0. The space of 
solutions 5? to the Klein-Gordon admits a natural symplectic form defined by 



*A<i>i,<h) ■= J tfWIbl n"(0iV^ a - 2 V„0i) . (2.14) 
s 

The integral is evaluated on a Cauchy surface S, but o is independent of the particular 
choice of E. This is due to the conservation, V^j^ = 0, of the current := (f)S7^(f)' — (p'V 1 ^ 
as can be shown by the application of Stokes' theorem to a timelike cylindrical spacetime 
volume bounded by portions of Cauchy surfaces. The space ay) may be called the 
covariant phase space of the theory. 

Alternatively, one can regard the canonical phase space. This is the space 

^:={(«,pU,pe^ t) K)} (2.15) 

of initial values on Ej of the Klein-Gordon equation equipped with the symplectic form 

<7^(Fi, F 2 ) := - J d 3 x (uipz - u 2 pi) (2.16) 



for Fi := (ui,pi) G i = 1, 2, where pi := A/|f)|(n M V M Mi) is the canonical conjugate to the 
configuration variable U{. 

The relation between the distinct phase spaces can be summarized as follows: By the 
unique correspondence between the solutions to the field equation and the initial values 
on a given Cauchy surface £ t , the spaces and are isomorphic, i.e., there exists an 
isomorphism 1 t : ^ — > 5? . This isomorphism induces a symplectic map (Jj(= T^ay, where 
X t * denotes the pulled back isomorphism. Consequently, both phases spaces are equally well 
suited for quantization |TV99] . 



2.2 Quantization 

Let us introduce some basic notions regarding algebras [BW92J. Let 21 be an algebra over 
C with a map * : 21 — >• 21 such that for all A,Be$l and a, (3 G C we have (a A + (3B)* = 
aA* + ]3B*, (AB)* = B*A*, and (A*)* = A. Then, * is called an involution and 21 is called 
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an involutive algebra or a *-algebra (star-algebra). If 21 contains a unit element 1 such 
that tA = At for all A G 21 then it is called a unital *-algebra. If the *-algebra 21 is also 
a Banach space where the norm satisfies ||AB|| < ||A||||i?|| for all A, B G 21 then 21 is a 
Banach *-algebra. If, in addition, \\A\\ 2 = \\A*A\\ then 21 is a C*-algebra. A *-subalgebra 3 
of 21 is a *-ideal if AB, BA G 3 for all A G 21 and B G 1 A C*-algebra 21 is called simple 
if it contains no non-trivial, i.e., different from and 21, closed *-ideals. 

The observables of the quantized theory are represented by the self-adjoint elements 
of a suitable algebra, e.g., a (unital) *-algebra or, if a stronger structure is desired, a C*- 
algebra. So, owing to the different necessities, there exist several formulations of the algebra 
of observables, which are substantially equivalent, albeit technically inequivalent. In this 
section, we introduce the formulation in terms of a Weyl algebra, which is a C*-algebra. For 
a thorough treatment of the quantization scheme in curved spacetimes see [Wa l94l |KW91| . 
All relevant facts regarding Weyl algebras and their representations can be found in [BR97J. 

Based on each of the symplectic spaces defined in section I2.1.3| we can define an abstract 
C*-algebra that obeys the canonical commutation relations (CCR) via the Weyl construc- 
tion. Consider a real symplectic vector space (Y, a). 

Definition 2.2 A Weyl algebra 2U('^ / ,(7) is a simple C*-algebra with unit generated by 
objects W(f) that are labeled by functions / G "f and that satisfy the relations 

(i) W(0) = 1, 

(ii) w(fy = w(-f), 

(iii) W{h)W{f 2 ) = e-^f^Wif, + f 2 ) 
for all fJ 1 J 2 EV. 

Condition (iii) is the Weyl form of the canonical commutation relations (CCR); Thus, a Weyl 
algebra is often called a CCR algebra. The elements of 22J represent the basic observables of 
the quantum theory. They are bounded operators, which avoids possible domain problems, 
and they correspond, formally, to exponentiated field operators W(f) = e -1 *^. This 
interpretation is mathematically well defined in regular representations (see section [2.2. 2p . 
Note that, since "V is not given a topology, the elements W(f) need not be continuous. 
Provided that a(f,g) is non-degenerate, the Weyl algebra is unique in the sense that given 
two Weyl algebras 2Ui and 2U2 there exists a unique ^isomorphism a : 2Ui — > W2 such 
that for any Wi G 2Ui and W 2 G 2B 2 , we have a-Wi = W 2 - 

States 

A state uj on the algebra of observables 21 is a positive, u>(A*A) > , normalized, lu(±) = 1, 
linear functional u : 21 — » C for all A G 21. The set of all states is a convex set, i.e., a 
mixture to := \\0Ji + • ■ ■ + \ n uj n of states u%, . . . , u n with \ > 0, A« = 1 is again a state. 
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A pure state is extremal in this set in the sense that it cannot be expressed as the sum of 
two other states with positive coefficients Aj. 

The general definition of states on an algebra is, on the one hand, clear and concise, but 
on the other hand, far to general for concrete applications. The space of states satisfying 
these conditions is enormous and requires further criteria that single out subspaces of states 
that are appropriate for a given physical situation. On a Weyl algebra one typically restricts 
attention to the class of regular states, which allows to introduce the quantum fields <&(/) 
as self-adjoint generators of the Weyl elements W(f). A further condition is to require 
the states to be quasi-free, i.e., to be completely specified by their two-point distribution, 
which makes the set of states tractable without removing most of the physically interesting 
states. Further conditions that we will use are homogeneity and isotropy, and the Hadamard 
condition. All these requirements will be introduced in due place. 

The algebraic approach makes contact with the traditional Hilbert space formulation of 
quantum theory via the GNS theorem, which says that every state uj on a C*-algebra 21 
gives rise to a representation of 21 on some Hilbert space. 

Theorem 2.3 (GNS construction) Let uj be a state on a C*-algebra 21 with unit- 
element. Then, there exists a complex Hilbert space J^, a unit-preserving representation 
7i u in terms of linear operators on Jf? , and a vector Vt w e such that 



for all A G 21. The vector is cyclic, i.e., 7r u ($l)Q u is dense in M > u . The representation 
tt^ is unique up to unitary equivalence. 

The triple ir w , ^) is called the GNS triple and the representation is called the GNS 
representation. 

The folium of a state uj is the set of all states that can be represented as density matrices 
p in the GNS representation of uj. So the folium consist of all states of the form 



with a positive, trace class, i.e., Trp < oo, operator p. An important theorem due to 
Fell [Fel60] states that the folium of a faithful representation of a C*-algebra is weakly 
dense in the set of all states. Since every physical experiment consists of a finite number of 
measurements and, furthermore, these measurements have limited accuracy, it is impossible 
to determine more than a weak neighborhood of a state. Thus, by Fell's theorem, it is 
impossible to find out in which folium the state lies. Note that, since all Weyl algebras are 
simple |Sim72j . all their representations are faithful. For a discussion of further implications 
of Fell's theorem see |Haa96l IWal94j . 

The GNS representations of different states need not be unitarily equivalent. In fact, 
the Stone-von Neumann uniqueness theorem fails for systems with infinitely many degrees 
of freedom and it is known that infinitely many inequivalent, irreducible Hilbert space 



uj{A) = (n u ,ir u (A)Q u )^ 



(2.17) 



uj(A) = Trpn(A) . 
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representations of the Weyl algebra exist (see, e.g., [Wal94|). Consequently, the folium of 
a single state does not encompass all possible algebraic states. This is in contrast to the 
finite case, where all irreducible regular representations, in particular, the Schrodinger and 
the Heisenberg representation, are unitarily equivalent. 

A representation tt on of the Weyl algebra 2U is called regular if the unitary groups 
A i — ► n(W(\f)), A G R are strongly continuous for all /. If tt is regular, one can, by 
Stone's theorem, introduce self-adjoint infinitesimal generators $n(f) of the Weyl elements, 
which act on Jif n . These operators may then be used to define annihilation and creation 
operators (see theorem 12.51) . A state lj on the Weyl algebra 2U is said to be regular if its 
GNS representation is regular. 

An automorphism a on a *-algebra 21 is a one-to-one linear mapping of the algebra onto 
itself that is compatible with the algebraic structure, i.e., it satisfies a(A ■ B) = a (A) ■ a(B) 
and a(A*) = a(A)* for all A, B G 21. A classical symplectic transformation on a) is 
a map that leaves the symplectic form invariant, i.e., a symplectic transformation is given 



by an operator T : Y — > "V such that a(T fi,T f 2 ) = cr(fi, fi) for all /i,/2 G Y. A 



symplectic transformation on a classical symplectic vector space corresponds directly to 
an automorphism on the associated Weyl algebra. For example, the time translation on 
stationary spacetime is implemented on the classical phase (^,o>) by a one-parameter 
group of symplectic transformations % : 5? — > =5^, which gives rise to a one-parameter 
group of automorphisms 



for all <p G 5? . We note that a pair (21, a t ) of a C*-algebra 21 and a strongly continuous 
automorphism group {a t } tgR acting on 21 is called a C*-dynamical system. This kind of 
system provides the basis for the definition of KMS states (see section 15.11) and passive 
states (section f4.1.ip . 

Given a Hilbert space representation tt of 21 on some Hilbert space Jif, we say that the 
symplectic transformation T is unitarily implementable if there exists a unitary transfor- 
mation U : Jif — > Jif such that 



for all A G 21. While there is no problem with the unitary implementation of time- 
translations in stationary spacetimes, the situation changes significantly for non-stationary 
spacetimes. The two-parameter family of symplectic transformations T t2ttl , describing time- 
evolution in that case, gives rise to a family of automorphisms a t2jtl on the algebra. However, 
these automorphisms are no longer implementable as unitary operators on a Fock space, 
as it has been shown in |TV99] for the Klein- Gordon field on the torus T 2 with non-flat 
Cauchy surfaces. One may actually conjecture that only transformations defined by the 
isometry group of a spacetime can be represented as unitary transformations on a Hilbert 
space. 



a t :W^W, a t (W ((/>)) :=W (%(/>)) 



(2.19) 




(2.20) 
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2.2.1 Hilbert space quantization 

In this section, we review the quantization of a linear, scalar field in a formalism that 
is close to the traditional Hilbert space quantization and directly applicable to curved 
spacetimes. The formalism starts with a real, symplectic vector space - in our case, the 
vector space (J?, ay) of solutions to the Klein-Gordon equation. Then, on ay) an 
inner product /i with suitable properties is chosen, which gives rise to a map K from 5? 
to a (complex) Hilbert space ffl . The quantum field theory is then constructed on the 
symmetric Fock-space J^^ 5 ) associated to the one-particle space M' '. We do not give all 
details of the construction, just the general procedure. The authoritative references caring 
for all contingencies are [KW911 IWal94j . 

First, we need to construct an inner product structure on the real, symplectic vector 
space(^, ay). So, choose any positive, symmetric, bilinear map \i : 5? x 5? — > R such that 

\o^u fo? < K<h, 0i)M02, <h) (2-21) 

for all 0i,02 6 . Since ay is non-degenerate, the map fi defines a real inner product on 
the vector space 5? . 

One can show that there always exists a fi satisfying ( 12.211) . but in a general curved 
spacetime there is no way to select a preferred one. While in the case of theories with 
finitely many degrees of freedom different choices of \x lead to unitarily equivalent theories, 
in the non-finite case the theories turn out to be, in general, unitarily inequivalent. In a 
stationary spacetime an operator K associated to \i can be defined that projects solutions 
in y onto the subspace of 'positive frequency solutions'. These solutions have positive 
frequency in a generalized sense, namely, with respect to the timelike Killing vector field £^ 
present in a stationary spacetime [AM75, Kay78, Wal94]. 



For non- stationary spacetimes there is no unique subspace of positive frequency solutions 
on which K could project. Nevertheless, one can proceed with a non-unique decomposition 
of y by the following results due to Kay and Wald. 

Theorem 2.4 ([KW91]) Let y be a real vector space on which are defined both a bilinear 
symplectic form ay and a bilinear positive symmetric form /i satisfying ( I2.2ip . Then, one 
can always find a complex Hilbert space J$? together with a real-linear map K : 5? — > J%? 
such that 

(i) the complexified range of K , i.e., K, c / — iKy, is dense in Jfc ', 
(li) yu(>i, (j) 2 ) = 3? {Kfa, Kfc)^ for all fa, (fr 2 G 5?, 
(iii) ay{fa, h) = 23 (Kfa, K<p 2 )^ for all fa, fa e 5?. 

Moreover, the pair (K, is uniquely determined up to equivalence, where we say (K', Jff") 
is equivalent to (K, Jrf?) if there exists an isomorphism V : Jff —>■ 3tf" such that VK = K' . 
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So, to every triple (y, <iy, /i) there corresponds a pair (J^, K). The pair (^, K) is called 
a one-particle Hilbert space structure. The equations given in theorem 12.41 are often stated 
in the form 

(Kfa, Kfa)^ = ftfa, fa) + ^ay(<Pu fa) • (2.22) 

A corresponding operator K : 5? —>■ Jif can be defined which projects into the subspace of 
'negative frequency solutions', where J%? is the complex conjugate Hilbert space to JF. To 
remind the reader, the complex conjugate space Jrf? differs from ffl by the scalar multipli- 
cation: c© / = c ■ /, / G J^, c G C, where the bar denotes complex conjugation. One may 
as well say that an antilinear isometry V satisfying T 2 = 1 makes the transition between 
the spaces Jif and J$?. It follows immediately that K + K = 1. 

Once having defined a one-particle structure K, the remaining quantization procedure is 
straightforward. Define by 

#" S (JT) :=C© JT©(jr® s (2.23) 

the symmetric Fock space over the one particle Hilbert space J^. To each solution in y a 
corresponding operator <Ty($>, •) on the Fock space ^ S {J^) is defined by 

(f)) ■= ia{K(j)) - ia*(K0) , (2.24) 

where the standard creation and annihilation operators a*, a on ^ s (Jt?) satisfy 

[a(V0,a*(VO] = (^%r > ( 2 -25) 
a(V>)fl = (2.26) 

for all if), if)' G M' and the vacuum state Q := 1©0©0©. . . . This concludes the construction 
of the quantum field theory. 

The quantization scheme presented here, leads directly to a Fock space representa- 
tion of the algebra. However, as already noted, these representations need not be uni- 
tarily equivalent. If /ii ^ /i 2 the resulting representations {JP S (J$?[), a^ x {^, </>)} and 
{^5(^2), cr,y^ 2 ($, 0)} may be unitarily inequivalent. It follows immediately that different 
notions of 'particles' arise by different choices of /1. In a spacetime without time translation 
symmetry no preferred choice for \x exists. On stationary spacetimes, owing to the existence 
of a timelike Killing field, a satisfactory definition of a preferred /i, and thus a meaningful 
notion of 'particles' can be given |AM75l Kay78| . 



The viewpoint taken in algebraic quantum field theory is that unitary equivalence on the 
level of representations is not fundamental to the quantum theory. Rather, the algebraic 
relations between the collection of operators {ay($,(f))} on ^ s (Jif) are important. Hence, 
one postulates that a net of local algebras satisfying these relations completely determines 
the quantum theory and the fields play a role similar to coordinates in differential geometry: 
useful tools for daily work, but dispensable for the essential assertions. 
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Before we come to the algebraic formulation, let us note that there is an alternative 
procedure for the selection of a representation. Rather than specifying the bilinear form 
fi, one may define a complex structure J, i.e., a bounded linear map satisfying J 2 = — 1, 
on 5? for which — cr(0, J<p) is a positive definite inner product. This procedure is, e.g., 
used in the geometric quantization programme |Woo 92j. Its application and importance 
for quantum field theory in curved spacetimes is illustrated best in a series of papers by 
Ashtekar and Magnon- Ashtekar [AM751 IA"MA"80al IA"MA8 0bj. Both schemes, the choice of 
an inner product /i and the choice of a complex structure J are roughly equivalent. See, 
once more, [Wa l94j for details. 

2.2.2 The algebra of observables 
The Weyl algebra 

We have seen how a Fock space representation {JP S (J^), o>>($, </>)} of the algebra of ob- 
servables can be obtained via the choice of a bilinear form fi(<ft,<j)). Now, a Weyl algebra 
(see definition I2.2p can be defined through the unitary operators 



The C*-completion, which is known to exist, of the space generated by all W(<p) via formal 
finite sums £V \W{(j)i) comprises a Weyl algebra, which can be seen as the minimal algebra 
of observables of a quantum field theory in curved spacetime. This makes sense, since 
the Weyl algebras arising from different choices of fi are always isomorphic even if the 
corresponding Fock space representations are unitarily inequivalent |Sla72j . 

A second Weyl algebra, which is isomorphic to the former via lemma 12.11 can be con- 
structed as follows. Define using the third property of lemma l27fl for each / e a smeared 
operator on J^ S (J^) by 



W(<j>) := e~ ia ^>® . 



(2.27) 



These operators satisfy the Weyl relations 

w{<f>y = w(-cf>) , 

W(0i)W(0 2 ) = 6-2^(^)^(0! + <j> 2 ) . 



(2.28) 
(2.29) 




(2.30) 



W(f) := e"^ , 

where we set W(f') := W(f) if Ef = Ef again define a Weyl algebra, 

w(fy = w(-f) , 

W{h)W{f 2 ) = e-^^Wih + f 2 ) , 



(2.31) 



(2.32) 
(2.33) 



where 




(2.34) 
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The field algebra An abstract definition of a *-algebra 21 of field operators for a real, 
linear, scalar field can be obtained as follows. Take the unit 1 and the formal smeared field 
operators <&(/), where / G and demand that 

(i) $(/) is linear, 

(ii) $(/) is hermitian: $(/) = $(/)*, 

(hi) $(/) satisfies the field equation: <&(p 8 + m 2 ]/) = 0, 

(iv) $(/) satisfies the commutation relations. [$(/i), $(/ 2 )] = iE(fi, f 2 )l. 

To be precise, one takes the free algebra over the field of complex numbers generated by the 
symbols $(/), $(/)*, and 1 and divides by the *-ideals generated by the properties stated 
above. 

The formal correspondence between the field algebra and the Weyl algebra becomes 
mathematically well defined in regular GNS representations. For those representations, as 
already noted, the map A h- > n(W(\f) = 7r(e _1 *^ A ^), A G R defines a strongly continuous 
unitary group for every /. Thus, by Stone's theorem, the operators $(/) are self-adjoint 
generators of these groups. These generators satisfy the requirements for a *-algebra of 
fields by lemma 12.11 The collection of all such operators may, equally well as their Weyl 
counterparts, be interpreted as the collection of fundamental observables of the theory. The 
same is valid for the collection of <Jy($>, <p) in a regular representation. 

2.2.3 Quasi-free states 

There is a distinct class of regular algebraic states, namely, the quasi-free states, also known 
as Gaussian states, which have GNS Hilbert spaces that look like the familiar Fock spaces 
built on a one-particle Hilbert space (see theorem 12.51 below) . The class of quasi-free states 
contains the usual vacuum states in stationary spacetimes as well as other vacua obtained 
by mode decomposition of the field operators. In quantum statistical mechanics, quasi- 
free states represent the general form of equilibrium states for free bosonic systems |BR97t 
HR97j. Besides this, quasi-free states are well suited for calculations as they are exclusively 
determined by their two-point distribution. 

A quasi-free state can be defined as an abstract linear functional of a special kind on 
the algebra of observables, or in a regular GNS representation by the requirement that it is 
completely fixed by the two-point distribution. Let us investigate the first possibility. To 
define a state on a Weyl algebra 2U, it suffices to specify its expectation values on the Weyl 
operators W G 2U. As before, let \i be a real scalar product on and consider the Weyl 
algebra 2U(J^, cry). Define a functional uj^ : 2B — > C by 



(2.35) 



2.2. Quantization 



25 



for all <p G 5? . Its action is extended to the whole algebra by linearity and continuity. Now, 
if \i satisfies ( I2.2ip then uj^ is positive on the whole algebra and thus a state. Any quasi-free 
state uj^ can be realized in a representation as the vacuum state in a Fock space by the 
following theorem. 

Theorem 2.5 ([KW91J) Let (K,J^) be the one-particle Hilbert space structure obtained 
fromuj^ by Theorem \2.4\ The GNS-triple (^C> M , 7Tu>uj ^> m ) of the state is equivalent to a 
triple (JF(J4f), tc^ , Q, ) with the following properties. 

(i) The GNS space ^(Jif) := J^ s (Jf?) is the symmetric Fock space built on J$f. 



(ii) The representation n 1 ^ is specified by ir^(W((f))) = e ^ K ^> a *( K <t>)] ; where the bar 
denotes the closure. 

(Hi) The state := 1 © © © . . . is the (cyclic) Fock vacuum in ^ S {J^). 

The purity of uj^ is equivalent to the irreducibility of the representation 7r . Moreover this 
is equivalent to the property that K<5? alone, rather than Ky + iKy, is dense in J$? 
|KW91j . 

The two-point distribution of a quantum field in the state u> on the Weyl algebra 
%3{y, ay) is defined by 



(2.36) 

s=t=0 



provided that the right hand side exists. For quasi-free states, as defined by (12.35)) the 
two-point distribution always exists and is given by 

uWtyufa) =y.(<t>iM + \vA<h,<h) • (2-37) 

If one works with the field algebra of operators $(/), a quasi-free state u can be defined, by 
the requirement that all odd n-point distributions vanish and the even ones are determined 
by the two-point distribution via the the recursion formula 



*)(d>(A), • • • , Hhj)) := EII w(2) (/4)./^)) (2-38) 

a i=l 

for j G N, where the sum is taken over all permutations a of {1,2, . . . ,2j} with cr(l) < 
er(2) < • ■ • < a(j) and a(i) < a(i + j),i = 1, . . . , j. 

According to theorem I2.5j the two-point distribution S(fi,f 2 ) of a quasi-free state can 
be calculated in the GNS representation by 

S(f lt f 2 ) := (^^(frWm*) (2.39) 
= (KEf 1 ,KEf 2 ) je (2.40) 

= t i(Ef 1 ,Ef 2 ) + ^a(Ef 1 ,Ef 2 ) (2.41) 
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Thus S(f, f) = n(Ef, Ef) = (KEf, KEf) ^, for / e 3>(M) and we see that we can define 
the generator (12.351) of a quasi-free state either by the use of /i, S, or ( , ) M> . 

Saturating the inequality (I2.2ip turns out to be equivalent to the irreducibility of the GNS 
representation arising from uj^. Thus, (remember the comments after theorem I2.5D such a 
\i results in pure quasi-free states. However, already in stationary spacetimes, thermal 
equilibrium states at finite temperature are represented by quasi-free states which fail to 
satisfy the saturated version of (I2.2ip . and in a spacetime with a non-compact Cauchy 
surface, states in the folium of a thermal equilibrium state do not, in general, lie in the 
folium of any pure quasi-free state |Wal94j . Thus, in order to incorporate thermal states, 
we cannot restrict ourselves to the case of pure states. Note that the common notion of a 
'vacuum state' or, more generally, ground state corresponds to a pure, quasi-free state in 
the algebraic approach. We note that all pure, quasi-free state are Fock states and thus 
related through a Bogoliubov transformation [MV68J. 



2.2.4 Hadamard states 

Even in the restricted class of quasi-free states there exist many states that cannot be con- 
sidered physical; this should not be misunderstood as the statement that all physical states 
are quasi-free. A further condition that is believed to reasonably narrow down the class of 
states is the Hadamard condition. It has several formulations and was used already for a 
long time (see, e.g., [DB60j ) before Kay and Wald |KW91j put it on a sound mathematical 
foundation. The Hadamard condition essentially restricts the singularity structure of the 
two-point distribution at coinciding points such that it comes close the Minkowski vacuum 
state. On a heuristic level, one could say that in a Hadamard state the high frequency 
modes of the field are 'close' to their ground state. 

The principle motivation for the Hadamard condition comes from the point-splitting 
renormalization scheme. The absence of a preferred vacuum state makes the normal ordering 
procedure of standard quantum field theory nonviable in a general curved spacetime. The 
point-splitting prescription is a replacement for normal ordering in the sense that it defines 
sensibly differences of stress-energy expectation values even in curved spacetime. 

The basic idea is as follows. The calculation of the stress-energy tensor involves the 
calculation of objects like (<f)(x) 2 ). These objects are, in general, ill defined as they involve 
the calculation of products of two distributions at a point. Such products are only well 
defined in special cases, where the wavefront sets (see below) of the involved distributions 
'fit' to each other |Hor03j . The point-split object {4>(x)(j)(x')) , however, makes sense as a bi- 
distribution on M.. Now, one demands that for physically reasonable states the singularity 
structure of ((f>(x)(f>(x')) as x' approaches x should be the same as for (O, 4>(x)(p(x')Q) . This 
provided, the difference (<fi(x)(j)(x')) — (Q,<f)(x)<j)(x')Q) is a smooth function, which allows 
taking the coincidence limit x' — > x. 
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States satisfying the Hadamard condition yield a renormalized stress-energy tensor T^ v 
that satisfies Wald's axioms |Wal94j . which are believed to be reasonable physical assump- 
tions on T^ v . In brief, the assumptions are as follows. For states u,Ui,U2- 

(i) (T^) — (T^) should be well defined by the point splitting procedure. 

(ii) T^v should be locally covariant. 

(iii) (T^) should be conserved, (T^) = 0, for all states uj. 

(iv) In Minkowski spacetime, (f2, T M!/ fi) = 0. 

These assumptions fix the expectation value (T^) uniquely up to a conserved local cur- 
vature term that is independent of uj. Actually, it is possible to locally construct a bi- 
distribution H(x,x') such that if 



is a smooth function for a state uj then (T M „) is well defined. We refer the reader to [Wal94j 
for subtleties of the construction. 

It has not only been proved that there exist many Hadamard states for the linear, scalar 
field in a globally hyperbolic spacetime; One also knows that the canonical ground state 
and the thermal equilibrium states in stationary, globally hyperbolic spacetimes are quasi- 
free Hadamard states |SV00] . Regarding interacting fields, the construction of renormalized 
perturbative quantum field theory in a generic spacetime is possible for Hadamard states 
[BF00, HWOT]. On the other hand, only few Hadamard states on curved spacetimes have 
been explicitly constructed. Among these are the ground states on ultrastatic spacetimes 
[FN W811 IJun96] and KMS states on ultrastatic spacetimes with compact Cauchy surfaces; 
Furthermore, we mention the adiabatic vacuum states of infinite order [Jun96lDun02j . None 
of these examples deals with thermal states on a non-stationary spacetime. 

The formulation by Kay and Wald 

Let t(x) be a any global time function that increases towards the future and let t(x,x') := 
t(x) — t(x'). Furthermore, let a(x,x') be the squared geodesic distance, i.e., 



where 7 is the unique geodesic connecting x and x' with parametrization y(-) and the sign 
is plus for spacelike y(-) and minus for timelike y(-). The geodesic distance is well defined 
and smooth in the set O in j\A x A4, which is taken to be a neighborhood of the set of 
causally related points (x, y) such that J + (x) fl J~(y) and J + {y) H J~(x) are contained 



F u> (x,x t ) := (<f>{x)(f){x')) ^ - H{x,x') 



(2.42) 




(2.43) 
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within a convex normal neighborhood. A convex normal neighborhood of a point x is a 
neighborhood U of x such that there exists a unique geodesic connecting each pair of points 
in U and staying entirely within U. A manifold can always be covered by such neighborhoods 
[HE75] . 

For each n G N and e G R, e > define a complex-valued function in (9 by 

= 7^ f^^r+t;W(x I x')lnaf(x,x')) , (2.44) 
(2tt) 2 V "*^^ ) / 

where 

<r*(z, x') = tr(z, x') + 2iet(x, x') + e 2 . (2.45) 

The smooth function u(x,x') is the square root of the van Vleck-Morette determinant, and 
v^ n '(x, x') G £(0) is a real function defined by the power series 

n 

t> (n) (x, x) := v ™{x, x>(x, x') m , (2.46) 

m=0 

where the determined by the Hadamard recursion relations. The branch cut 

for the logarithm is taken to lie along the negative real axis. The field equations and 
commutation relations require that u and v are uniquely determined by the local geometry, 
i.e., by the metric and its derivatives. Of course, G t e ' n (x, x') is singular for coinciding 
points x and x', but it is also singular for points connected by a null geodesic. 

Let E be a Cauchy surface and let M be a causal normal neighborhood of S [KW91J. 
Denote by O' an open neighborhood in jV x J\f of the set of causally related points such that 
the closure of O' in J\f x M is contained in O. Now, define a function x{ x i x ') ${N x A/") 
with the property 

for (x, x ) G O . 



A state a; is said to be a Hadamard state if its two-point distribution satisfies the following 
requirement: for each n G N there exists a C n -function w n (x,x') on M x A/" such that for 
all /1J26 5^ (A") we have 

5(A,/ 2 ) = lim J dfi(x) d^{x')h{x)h{x') Af (x,x') , (2.48) 
AfxM 

where 

A*' n (x, x') = x(x, x')G*' n (x, x') + w n (x, x') (2.49) 

and the measures are the induced measures <i/i(x) := g^xa/Ijj]. The functions w n (x, x') 
depend on the individual state in the folium. Note that the function x(x,x') introduces a 
localization of the singular object G*' n (x, x') in spacetime, which has the consequence that 
the Hadamard condition only cares for the ultraviolet modes of a field. 
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The microlocal formulation 

For linear quantum fields the microlocal spectrum condition, which we will state below, 
is equivalent to the Hadamard condition. This has been shown by Radzikowski |Ra d92. 
IRad 96] . The microlocal spectrum condition has been investigated and extended to curved 
spacetimes and Wick powers of scalar fields in [BFK96J. In this sense, it is more general 
than the original Hadamard condition as it also can be generalized to non-linear fields. 
Microlocal analysis shifts the study of singularities of distributions from the base space to 
the cotangent bundle. The 'micro-localization' property makes this formulation well suited 
for curved spacetimes. 

We define the notion of wavefront sets by first introducing the set of regular directed 
points. Let u G @'(R n ) be a distribution. A point in the cotangent bundle of R n is 

called a regular directed point of u if there exists a smooth function <p G @(X), X C R n 
which does not vanish at x and such that for any m G N there exists a constant B m satisfying 

\Me)\<B m (i + \e\r m (2.50) 

for all £' in a conical neighborhood U C R n \{0} of £. A neighborhood U is called conical if 
£ G U t£ G U,t G R + . Now the wavefront set WF(u) of a distribution u G &{X) is the 
complement in X x R n \{0} of the space of all regular directed points of u. 

Theorem 2.6 (Microlocal spectrum condition) A quasi-free state u of the Klein- 
Gordon field on globally hyperbolic spacetime Ai is a (global) Hadamard state if and only if 
its two-point distribution has the wave-front set 

WF{J 2 \x 1 ,x 2 )) = C + , (2.51) 

where 

C + = {[(a*, £i),(s3, 6)] G (T*(M)\0) x (T*(M)\0);(x 1 ^ 1 ) ~ (x 2 ,-£ 2 );£? > 0} . (2.52) 

The notation (xi,£i) ~ (^2^2) means that x\ and x 2 can be connected by a null geodesic 
such that £f is tangential to 7 at X\ and £ 2 is the parallel transport of £f along 7 at x 2 - 

So, singularities in the two-point distribution lj( 2 \x, x') only occur if x and x 1 are lightlike 
connected and the singularities have positive frequencies. 

Wavefront sets have the property that for two distributions u, v it holds that 

WF(« + !))CWF(«)UWF(d). (2.53) 

Thus, WF(m) C WF(m -v)U WF(t») and WF(u) = WF(t») if WF(u - v) = 0, i.e., the 
distributions u and v have the same wavefront set, and hence the same singularity structure 
if the difference of u and v is a smooth function. 
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3 Scalar fields in Robertson-Walker 
spacetimes 



Robertson- Walker spacetimes are a particularly important class of spacetimes to investigate. 
They are homogeneous and isotropic solutions of the Einstein equations. These assumptions 
on the geometry of spacetime determine the solutions up to three discrete types of spatial 
geometry - flat, closed, and open geometry - and an arbitrary positive function a(t), which 
determines the expansion of the spacelike sections. Although these spacetimes do not 
possess a time evolution symmetry, they possess a preferred class of observers, namely, the 
ones who are orthogonal to the homogeneous surfaces. These are also called isotropic (or 
comoving) observers. Owing to the present symmetries, the field can be written, in the usual 
way, as an integral (or sum, for the closed geometry) over modes. The spatial dependence of 
the modes is known explicitly, but the time dependent part does not acquire the usual form 
e ±tuJt . Rather, it satisfies a certain differential equation with time-dependent coefficients. 
A determination of solutions to this equation is possible by a method motivated by a 
WKB type approximation. This led Parker to define the adiabatic vacuum states [Par69j. 
The adiabatic vacuum states are defined such that the particle creation is minimized in 
an expanding universe. Later, the adiabatic vacuum states were redefined by Liiders and 
Roberts [LR9 0] in a more rigorous setting of quantum field theory in curved spacetime. 

In this chapter we introduce, first, the quantum theory of the scalar field in the formu- 
lation used in |LR90j . We also give a brief account of adiabatic vacuum states. Then, we 
quote the important theorem on the general form of the two-point distribution of a ho- 
mogeneous, isotropic, quasi-free state for the scalar field in Robertson- Walker spacetime. 
Finally, we calculate, as a preparation for the construction the almost equilibrium states, a 
'four-smeared' version of the two-point distribution. 

3.1 The algebra and the states 

We have seen in chapter [2] how the Weyl algebra or the field algebra for a quantum field 
is constructed in a general, curved spacetime. There is an alternative way to obtain the 
algebra of observables, which goes back to Araki, namely, by constructing a self-dual algebra 
|Ara68t IAS72} IAY82j . The self-dual algebra is a different route to the construction of the 
CCR algebra, distinguished by the fact that one chooses a different set of generators that 
act on a complex vector space JtT. This approach is used in [LR90j , where the Weyl algebra 
associated to a self-dual algebra is is taken as the algebra of observables. 
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The self-dual algebra 

A self-dual algebra is based on a phase space triple (J?T, T, 7), where is a complex linear 
space, T is an antilinear involution of J?T, and 7 is a Hermitian form on satisfying 

1 (Tf,Tg) = -~ r (f,gy . (3.1) 

The space Jt^ serves as the test function space for creation and annihilation operators. 
The indefinite inner product ~y(f, g) gives rise to the canonical commutation relations and 
the involution T combines complex conjugation with the interchange of test functions for 
creation and annihilation operators [AY82]. 

For the Klein-Gordon field on a globally hyperbolic spacetime (Ai,Q^ u ) one defines 

JT := &{M)I[{U S + m 2 )£l(M)] , (3.2) 
j(f,g):=iG(f,g) , (3.3) 

where G(f,g) := J d 4 x fEg, f,g& J^, and T is defined by the requirement 

17 = / ■ (3.4) 

Note that @(M) := C^(Ai) is a space of complex functions, as opposed to the space of 
real functions used in the construction in chapter [2j and the fundamental solution E is 
defined in equation (12. 7L Equation ( 13 . 21) ensures that the corresponding field satisfies the 
Klein-Gordon equation. The involution V satisfies (13.11) by the properties of the Green's 
function G(f, g). 

On the phase space ( J^, T, 7) one constructs the self-dual CCR algebra 2l( J?^, T, 7) as 
follows. First, consider the complex, free, *-algebra over J£T generated by the symbols 4>(f), 
their conjugates 0(/)*, and the identity 1, where / G J^. Then, take the quotient of it by 
the two-sided *-ideal that is generated by the relations 

4>{f) is complex linear in / , (3.5) 

W)>(9)-^(/r = 7(/,9)l, (3.6) 
<P(Tfy = 4>{f) . (3.7) 

for f,g G In order to circumvent problems with unbounded operators, one passes, as 
usual, from 2l(J?T,r,7) to an associated Weyl algebra. The Weyl algebra 22J(J^R,i7) is 
based on the real linear space :— {f G :Tf — /} equipped with the real symplectic 
form vy(f,g): 

W(f)W(g) = e-WMw(f + g) . (3.8) 

Note that this is formally equal to (12.321) if J^sr coincides with £T . A net of C*-algebras is 
obtained by taking %{0) to be the C*-subalgebra generated by the Weyl elements W(f), 
f G @(M) with supp / C O. 
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Quasi-free states 

A quasi-free state for a real, linear, scalar field is specified by its two-point distribution 

S(f,g)=uj s (<J ) (fy<j ) (g)) (3.9) 

for f,g G J£T. This two-point distribution may be seen as a complex scalar product on the 
space . It specifies the corresponding quasi-free state in terms of Weyl operators by 

u S (W(f)) := e-**C/./> . (3.10) 

The two-point distribution of a quasi-free state is a polarization on the phase space 
T, 7). This means that S(-, ■) is a positive Hermitian form on J£T such that 

S(f,g)-S(Tg,Tf) = ~r(f,g) . (3.11) 

Given S(-, •) one defines a scalar product (•, -)s on by 

(f,g)s = S(f,g)+S(Tg,rf) (3.12) 

so that by the Schwarz inequality and the triangle inequality one obtains |AS72] 

\i(f,g)\ 2 <(fj)s(g,g)s, (3.13) 

which is the complex version of (12.211) . 

It is remarkably easy to characterize pure, quasi-free states in this scheme. Namely, 
denote by Jffs the Hilbert space completion of J?T by (-,-)s- On we have an operator 
5" satisfying (/, Sg)s = S(f, g). The state us is a pure quasi-free state if and only if S is a 
basis projection on J(fs |AS72] . 

3.2 The spacetime and the field 

As already noted, Robertson- Walker spacetimes are homogeneous and isotropic solutions to 
the Einstein equations. They are globally hyperbolic Lorentzian manifolds with topology 
Ai £ — R x S £ , where e = 0, +1, —1 discriminates three types of spacelike hypersurfaces. The 
Cauchy surfaces S e are homogeneous Riemannian manifolds with constant curvature of sign 
e. The homogeneous and isotropic spacetimes can be endowed with the Robertson- Walker 
metrics 



ds 2 = dt 2 - a(tf 



-^--+r 2 (d6 2 + sm 2 e d V 2 ) 
1 — er 2 v ' 



(3.14) 



where the coordinates cover the ranges r 6 [0, oo) , 6 G [0, 7r] , (p G [0, 2n] for e = 0, — 1 and 
r G [0, 1] for e = +1. The function a(t) is a strictly positive, smooth function describing 
the expansion of the universe, and H(t) = ^ is the Hubble parameter. 
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The induced metric on the Cauchy surfaces can be written as 

ds 2 = dt 2 - a(t) 2 5 £ ij dx i dx j , (3.15) 

where we wrote the induced metric on the Cauchy surface as r)f„-(t) = a(t) 2 sfj. Note that the 
metric rjy is time-dependent while Sij is not. We use E to denote the manifold with metric 
Sij, while E 4 is endowed with the metric fjy. The future directed normal fields of all the 
hypersurfaces E e are given by n M = (1, 0, 0, 0). These fields are geodesic, i.e., n fl V^n l/ = 0. 
It is convenient to regard the Cauchy surfaces E £ as being embedded in R 4 by 

E° 
E+ 



Generally speaking, one calls the spacetime M. + a closed universe since E + is compact. It 
is also customary to call Ai~ and A4° open and flat universes, respectively. Each of the 
Cauchy surfaces E e is a homogeneous surface for a different symmetry group G e . In detail, 
these are the groups: G° = E(3), the Euclidean group; G + = SO(4), the rotation Group; 
and G~ = £+(4), the Lorentz group. 



The field equation 

In Robertson- Walker spacetimes the Klein-Gordon equation ( 12.2ft becomes 

2 

-^<f>{t, x) + 3H(t) -^<f>{t, x) + (-A„ + m 2 )0(t, x) = , (3.19) 

where Af, denotes the Laplacian on Ef. If one assumes that the field modes u^x) = 
Tk(t)lk(x) and their complex conjugates form a complete orthonormal basis, the general 
solution of the Klein-Gordon equation can be written as 



x G R 4 : x° = 

G R 4 : (x ) 2 + ^(x 1 ) 2 = 1 



x 



i=l 
3 



x G R 4 : (x ) 2 - ^(a;*) 2 = 1 ,x°>0 



(3.16) 
(3.17) 
(3.18) 



i=l 



0(t,x)= / dk[T k (t)y k (x)a(k)+T k (t)y k (x)a(k)*] 



(3.20) 



3.2. The spacetime and the field 



35 



We think of a(k), a(k)* as arbitrary complex coefficients with no reference to particles. This 
decomposition of the field is valid on all types of Robert son- Walker spacetimes if one defines 
the measure dk accordingly. Our notation is as follows: 

(h,k 2 ,k 3 ) G R 3 , k := |k| , (3.21a) 

(k,Z,m) G N x N x Z , (3.21b) 
(h,k 2 ,k 3 ) G R 3 , k := |k| . (3.21c) 

Note the subtle difference in notation between the absolute value of the three-momentum, 
k := |k|, and the four- momentum k := (ko, k\, k 2 , k^) = (fco, k). The functions T^(t) depend 
on k and the functions lk(x) depend on k. 



£ = : 

e = +1 : 
e = -1 : 



dk := / d 3 k , 

dk := 
dk := 



R 3 

oo k I 

EEE- 

k=0 1=0 m=-l 

d 3 k , 

R 3 



k := 

k := 
k := 



Homogeneous, isotropic, quasi-free states 

The groups G £ , act as isometry groups on the manifolds Ai £ := IR x S £ by g(t,x) = (t,gx), 
g G G e . It follows that they must commute with the fundamental solution E and thus act as 
a group of transformations on the phase space. In turn, this defines a group of (Bogoliubov) 
automorphisms {a g } on the Weyl algebra via 

a g (W(f)) := W(gf) (3.22) 

for g G G £ . Now, we say that a state u is homogeneous and isotropic if u o a g = u, 
g G G £ . Since a quasi- free state u s is, by definition, uniquely determined by its two-point 
distribution S(fi, f 2 ), this translates to the necessary and sufficient condition 

S(g/i,g/ 2 ) = S(/i,/ 2 ) (3.23) 

for all g G G £ . The G-invariance of the two-point distribution is analyzed in [LR90] by 
passing to the phase space (^#, a) of initial data for <j){x) at time t, where 

J: := { (u, a.(t) 3 p) , it, p G ^(E, R)} (3.24) 

and the symplectic form is a variant of ( 12. 16[) : 

&{F U F 2 ) := a(t) 3 / dx ( Pl u 2 - Ml p 2 ) (3.25) 
s 

for Fj := (uj,a(t) 3 pj) G i = 1,2, where dx = g? 3 :ea/Ts]. Please note that a is defined 
using the measure derived from Sij and not from t)^. 
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By introducing suitable topologies on one finds that it is sufficient to compute the 
commutants of the unitary representations U e of G £ on L 2 (E) given by Z7|/ :— f o g -1 , 
/ G i^(E) [LR90] . The representations U £ are decomposed into irreducible representations 
by (generalized analogues to) Fourier transforms 

/(k) := J dx F k (x)/(x) (3.26) 

for / G L 2 (E). In each case, the Fourier transform is a unitary operator from L 2 (E e ) to 
L 2 (E e ), where E £ denotes the momentum space associated to E e . And, again in each case, a 
bounded operator on L 2 (E £ ) commuting with U e corresponds on L 2 (E £ ) to a multiplication 
by a bounded function of k. 

The functions VL( X ) constitute an orthonormal basis of eigenfunctions of the Laplacian 
A s in L 2 (E) (see below for the precise forms). The inverse Fourier transform is given by 



/(x) = J dk r k (x)/(k) (3.27) 

and one has the usual completeness relations 

y dk Y k (x)F k (x') = 5(x, x') , (3.28) 
y dx F k (x)F k ,(x) = 5(k, kO • (3.29) 

The <5(k, k')-distribution is to be taken with respect to the measures dk defined in (13.21 j) : 

y dk7(k')5(k,k') = /(k) . (3.30) 



Spatial mode functions 

We give a brief account of the functions Y^ix.) in each of the three types of Robertson- 
Walker spaces. What is needed is a direct sum or direct integral decomposition of the space 
L 2 (E £ ) in terms of eigenfunctions of the Laplacian A s , i.e., 

A s F k (x) = -£(k)Y k (x) . (3.31) 

In each case, the decomposition is obtained by different means, but it exists and allows to 
treat, to a great extent, the solutions in the different types of Robert son- Walker spacetimes 
on an equal footing [LR90J. 
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[e = 0]: In the flat case we have Aj, = -^A s . Thus, it follows that A s Yk(x) = — k 2 Yk(x) 
and that the generalized eigenvectors lk(x) are independent of t. The eigenvectors and 
their negative eigenvalues are given by 

F k (x) = -^re ikx , E(k) : = k 2 . (3.32) 

The direct integral decomposition amounts to the ordinary Fourier transform. Note the 
useful relations lk(x) = Ylk(x) and |Yk(x)| : 



[e = +1]: For the closed universe, the solutions lk(x) are the spherical harmonics 

Y k (i/>,0,<p) = A kl n+(^) Y hm {#,<p) , E(k) :=k(k + 2) , (3.33) 

(k = 0, 1, . . . ; I — 0, 1, . . . , k; m — —I, — l+l, . . . , I), where the Y^ m are the harmonics on the 
2-sphere, the IlJ are real polynomials in sinip and cosip, and the A k i are real normalization 
constants. One has a direct sum decomposition L 2 (£ + ) = 0^1 o 7^k, where Ti k denotes the 
span of lk,z,m as I and m vary. 

[e = — 1]: In the open universe the Cauchy surface S~ is regarded as being embedded in 
M, £ = (1, £) G M and x-£ is the Minkowski scalar product. A set of generalized eigenvectors 
are 

^k(x) = - — — g-(a; ■ 0~ 1+ik , ^(k) := k 2 + 1 , (3.34) 



(27T) 



2 



where k = k£ G R 3 . Here, the Fourier transform is a map with values in L 2 (S 2 , dQ), where 
S 2 is the two-sphere embedded in R 3 , i.e., a function on the set of horospheres, in the 
language of jGGV66j . 

Time-dependence 

The time-dependent function T k (t), which appears in (I3.2(jp . is required to be a solution to 
the differential equation 

f k (t) + 3H(t)f k (t) + u>l(t)T k (t) = , (3.35) 
where the frequencies ui k (t) are given by 

w£(*):= ^ + m2 ( 3 - 36 ) 
and, additionally, to satisfy the condition 

T k (t)f k (t) - T k (t)T k (t) = -±- . (3.37) 
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The explicit form of E(k) in each type of Robertson- Walker spacetime has been given in 
the last section. The condition (I3.37p . which amounts to the Wronskian, determines the 
commutation relations of the a(k), a(k)* [HNS84IIPF74] . 

Equation (13.35!) i s an ordinary, homogeneous, linear differential equation with variable 
coefficients. In case of a static spacetime it has the explicit solutions 

T k (t) = -«^=e-^* . (3.38) 

In the general case it has a fundamental system T k (t),T k (t) of solutions, which cannot be 
calculated explicitly. 



Adiabatic vacuum states 

The only freedom one has in the construction of pure, homogeneous, isotropic, quasi-free 
states is the choice of initial data for the function T k (t). This choice can be made on 
physical grounds. For example, one can chose initial values such that the resulting states, 
in a certain sense, minimize particle creation in an expanding universe and reduce to the 
known particle notion in the static case. Then, one obtains the so-called adiabatic vacuum 
states |Par69j . 

In |LR90] the former definition of adiabatic vacuum states was put on a firm basis. The 
adiabatic vacuum states were redefined by fixing the large k behaviour of T k (t) and T k (t) 
such that the resulting states adhered the principle of local definiteness [HNS84J, which, 
roughly, requires that the the set of expectation values measurable in a bounded region O 
of the spacetime in a GNS representation of a state u should not depend on the state u. 
This rules out inequivalent representations of the local algebra of bounded observables 21(0) 
arising from different global states, which would be considered a pathology. For a precise 
definition and discussion of the principle of local definiteness and other related notions see 
|Ver94j . 

One can obtain iterative solutions T£(t) of (13.351) and hereby the adiabatic vacuum states 
|Par69j ILR90| . by using a WKB-type ansatz 

?k(t) = - 1/2 exp I i J dt n k (t) ] (3.39) 



(2a(tyn k (t)y 



to 



with yet undetermined positive functions £\(t). This ansatz satisfies the normalization, 
and it satisfies (I3.35P if 



2 



2 

!. V a. / 2 a. I \ fl k I 2 <\ 



^1 = ^1-- -) + 1 SO -~- (3-40) 



3.3. The two-point distribution 



39 



Now start with an iteration 

)(0) X 2 



«rj (3-4i) 



4 Va/ 2 a 4 I f)W / 2 o (n) 

The iteration procedure breaks down when we obtain a negative value for . This 

does not happen if one restricts to a time interval, t G /, where / G R is finite, and, in 
addition, k is chosen sufficiently large, k > k min |LR90j . Then, one defines for t ,t G / the 
function 

Wt\t) = -^exp [ i I dt n<»>(f) I . (3.43) 




An adiabatic vacuum state of order n is the pure, quasi-free state obtained as the exact 
solution of (I3.35P with the initial conditions 

T k (t) = Wl n \t) , f k (t) = Wl n \t) . (3.44) 



The adiabatic vacuum states depend on several quantities involved in their definition. First, 
they depend on the initial time t used for the initial values in (13. 44ft . This has no effect on 
the adiabatic vacuum state, as it amounts to common phase change of the initial conditions. 
Second, they depend on the extrapolation of £1^ to small k, which is always possible in 
a continuous manner, and amounts to some Bogoliubov transformation on the state, not 
affecting the large k behaviour. Of course, they depend also on the order of iteration, n. 

It has been shown in |Jun96l IJun02j that adiabatic vacuum states of infinite order are 
Hadamard states. This fact is of indirect importance for our construction as is is used 
in [Olb07a] to prove the Hadamard property of the states of low energy, and we base our 
proof of Hadamard property on the latter fact. The notion of adiabatic vacuum states was 
extended to more general spacetimes by the usage of Sobolev wavefront sets in [JS02] . 



3.3 The two-point distribution 

The two-point distribution of a homogeneous, isotropic, quasi-free state can be characterized 
by different means. For example, in [LR90J it is given in a form that takes data on the surface 
S. So, let Fi := (iii, a.(t) 3 pi) G ^ be a pair of initial values on the phase space associated 
to S. 
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Theorem 3.1 ([LR90]) The homogeneous and isotropic states on Robertson- Walker 
spacetimes have two-point distributions of the form 

J 2 \F,F>) := J dk(F(k),S(k)F'(k)) , (3.45) 

where 

i 

<F(k), S(k)F'(k)) := £ #<(k)fy(k)*}(k) . (3.46) 

i,j=0 

Here k i— > i%(k) is measurable and polynomially bounded. For the matrix S(k) it holds 
almost everywhere in k that 

5 i(k)-5 10 (k) = i (3.47a) 

S 01 (k) = S w (k) (3.47b) 

Soo(k) > (3.47c) 

Soo(k)Sn(k) > |S 01 (k)| 2 . (3.47d) 

Obviously equations (I3.47cl) implement (I3.47cp the positivity of the state (compare to 
(I2.2ip ). We remark that the two-point distributions depend only on the magnitude k of the 
three-momentum k because of the symmetry of the states. Exploiting these relations, it is 
possible to write 5*oi(k) = 5*^(k) + |. So, with J = (\ J) the matrix S(k) can be written 

as 

S(k) = S+^J, (3.48) 

where we introduced 

S := ( S JS f) ■ (3-49) 

We distinguish the original matrix 5'(k) and its real part S by omitting the k-dependence 
in the latter. This will cause no trouble since, from now on, we will almost exclusively deal 
with S. We remark that equation (13.47dl) becomes 

[S] > \ , (3.50) 

where [S] = det S. Equation (13 . 50 j) resembles a generalized uncertainty relation for the co- 
variance matrix of a state in quantum statistical mechanics. States of minimum uncertainty 
are pure states, which satisfy [S] — j. 

The two-point distribution of pure, quasi-free states is given by (13.451) with 
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where p(k) and q(k) are (essentially polynomially bounded and measurable) complex-valued 
functions satisfying 

W)p(k) - g(kMk) = -i . (3.52) 
The functions p(k) and g(k) are proportional to the initial data of T k (t) [Jim96j: 



S(k) 



a 6 TkT k -a 3 T k T k 
-a 3 T k T k T k T k 



The fundamental solution 



For the calculation of the explicit four-smeared two-point distribution of a quasi-free state 
we will need the following characterization of the fundamental solution E (see |LR90] ) . 
First, define the generalized function G(x,y) by 

G(x,y) := J dk G k (z >3 /°)y k (x)F k (y) , (3.53) 
G k (x°, y°) := i(T k (x°)T k (y°) - T k (x°)T k (y )) . (3.54) 



Then the fundamental solution E : 3s(M) — > $(M) can be written as 

(Ef)(x°,x) := -G(xJ) := J dy° J dk Y k (x)G k (x°, y°)f(y°, k) , 



(3.55) 



where 



/V, k) := a(y ) 3 | rf 3 y r k (y)/(y°, y) . (3.56) 

The kernel G k (-,y°) satisfies for each fixed y° equation (I3.35P (by linearity) with initial 
conditions 

G k (y°, y°) = i(T k (y°)T k ( y °) - T k ( y °)T k (y )) = , (3.57) 

G k (y°, y°) = i(f k (y°)T k (y ) - T k (y°)T k (y )) = . (3.58) 

This shows that G k (x°,y°) is independent of the particular solution T k (t) used in equation 
3.3.1 Four-smeared two-point distribution 

As a first step to our construction, we need the explicit form of the two-point distribution 
S(f, g) of a homogeneous, isotropic, quasi-free state in terms of the real matrix S from 
(13.491) and the solutions to equation (I3.35p . We prove the following lemma. 
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Lemma 3.2 Let lu be a quasi-free, homogeneous, isotropic state of the Klein-Gordon field 
in a Robertson- Walker spacetime M.. Then the two-point distribution of uj is given by 

S(f,g) = J dfi(x) J d^(x') f(x)g(x') J 2 \x,x') (3.59) 

with dfi(x) := v / 1flM 4a; > where 
uo {2 \x,x') := J dk F k (x)F k (x') (3.60) 
x T k (t)T k (t') U - + T k (t)T k (t') • b 2 + T k (t)T k {t') ■ b 2 + T k (t)T k (t') U + 

and 

h := S 00 \T k (t )\ 2 + &(t ) 6 S u \T k (t )\ 2 + (r k (t )f k (t ) + T k (t )T k (t )) a(t ) 3 ^ , (3.61a) 

b 2 := -S 00 Tl(t ) - a(to) 6 5 n tk(t ) - 2T k (t )t k (t )a(t ) 3 ^ 1 . (3.61b) 

Proof. In order obtain this formula for the two-point distribution, we need to calculate the 
two-point distribution of a homogeneous, isotropic, quasi-free state (13.451) . using the initial 
values F,H e ^# defined by 

a(to) 3 pi/) ' H := (Ah) (3 ' 62) 

on a Cauchy surface at time to- First, we need some Fourier transforms. Using equations 
(I3.26p . (13.291) and the abbreviated notation 

PoG k (x°,y ) :=G k {t ,y°) , (3.63) 
Pl G k (x ,y ):=G k (t Q ,y ) (3.64) 



one obtains 



Po£/(k) = J dy° G k (t , y°)f(y°, k) , (3.65) 

po^f(k) = J dy° G k (t , y°)!(y , y) , (3.66) 

^Ef(k) = J dy° G k (t ,y )f(y°,k) , (3.67) 

^Ef(k) = J dy° G k (t , y°)!(y°, y) . (3.68) 
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The Fourier transform of &(t ) 3 piEf is simply a(i ) 3 P\E f . Using this and (13.621) we can 
write the integrand of the two-point distribution (I3.45P as 

i 

F(k), S(k)tf (k)} = ^(k)^(k)^(k) (3.69) 

i,j=0 



p Ef(k)S 00 (k)p Eh(k) + a(t o y Po Ef(k)Soi(k) Pl Eh(k) (3.70) 



+ a (t ) 3 Pl Ef(k)S 10 (k)p Eh(k) + a (t ) 6 Pl Ef(k)Sn(k) Pl Eh(k) 
= J dp{y) J d»(y')f(y ,y)h(y' ,y')Y k (y)Y k (y')x (3.71) 

G k (t , y°)S 00 (k)G k (t , y'°) + a(t ) 3 G k (t , y°)S 01 (k)G k (t , y'°) 
+ a(t ) 3 4(to, y°)S w (k)G k (t , y'°) + a(t ) 6 G k (t , y°)S u (k)G k (t , y'°) 

where dp(y) := \/\g\d 4 y = a(y ) 3 ^/\s~\dy°d 3 y. We can now write the two-point distribution 



as 



S(f, h) = J dk (F(k), S(k)H(k)J (3.72) 

= J dp(y) J dp(y') f(y)h(y') uj^(y,y') (3.73) 

with 

^ 2 \y,y'): = J dk F k (y)F k (y') (3.74) 

x G k (to,y )G k (t ,y' )S 00 (k) + G k (to,y )G k (t ,y l0 )&(t fS 01 (k) 

+ G~ k (t , y°)G k (t , y'°) a (t fS w (k) + G~ k (t , y°)G k (t , y'°) a (t ) 6 S 11 (k) . 

We need to calculate the quantity in the square brackets. From the definition (I3.54p of 

G k (x°,y°) we have 

G k (t„, y°) ■= i (T k (t )T k (y°) - T k (t Q )T k (y )) . (3.75) 

In the following we use the abbreviated notation 

G := G k (t Q , y°) , G' := G k (t , y'°) , (3.76) 

G := d t0 G k (t , y°) , G' := d to G k (t , y'°) , (3.77) 

and 

T := T k (t ) , T := T k (y°) , T := T k (y'°) . (3.78) 
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Note that 

G = i(T T — T T) , 
G = i(T T — T T) , 

This gives the intermediate results 



Q = i(T V - TqT') , 
G' = i(T T 7 - f V) . 



GG 1 = \T \ 2 TV - t\tV - TqTT' + \T \ 2 TT' , 



• 2 



GG' = \T \ 2 TV - T Q TT 'obtain - T TT' + \T \ 2 TT' , 



GG' = T T TV - T T TV - T T TV + T T Q TT' , 



GG' = TqTqTT - T T TT' - T T TT + T T TT' . 
Making use of S i — + \ and Siq = — \ we finally arrive at 



[• 



= TV 
+ TT' 



S 00 \To\ 2 + a^nlTol 2 + (T T + T T )a 3 S, 



3 c» 
01 



— 'S'oo^o — a ^n-^o — 2T T a S' 01 



+ TV 
+ TT' 



—SqqTq — zPSuTq — 2ToToa 3 S , ( 



01 



S 00 \T \ 2 + a b S n |T | 2 + (T T + T T )aX + - 



(3.79) 
(3.80) 

(3.81) 
(3.82) 
(3.83) 
(3.84) 



(3.85) 



which gives the lemma. 



□ 



4 Quantum energy inequalities 



This chapter starts with a brief overview of quantum energy inequalities and their in- 
terpretation as stability conditions. Then, we quote the worldline quantum inequality of 
Fewster |Few00j . which we will use as the basis of our construction of almost equilibrium 
states. Afterwards, we give an explicit expression for the worldline averaged stress-energy 
in Robertson- Walker spacetimes in a form that is suitable for the minimization procedure 
to be accomplished in chapter [5j 



The geometry of spacetime is related to its matter content via Einstein's equation 



Without restrictions on the stress-energy tensor there would be no restrictions on the metric, 
thus no clue on physically realizable solutions. Apart from covariant conservation of energy, 
V M T^, the stress-energy tensor is believed to obey several more energy conditions, of which 
the most fundamental one is the weak energy condition (WEC). It says that 



for timelike is satisfied by all forms of physically reasonable matter. Interpreting w M as 
the four velocity of an observer, it guarantees that all observers at all points in spacetime 
always measure nonnegative energy density. 

The energy conditions are crucial ingredients in many important results concerning the 
behaviour of solutions to Einstein's equations. For instance, a result that requires en- 
ergy conditions is the positivity of the asymptotic gravitational mass of isolated objects 
|SY79[ ILV8H IWit81j . This has implications on the stability of Minkowski spacetime. Fur- 
thermore, energy conditions ensure that entertaining phenomena like traversable wormholes 
and stargates in 'designer spacetimes', time machines, i.e., spacetimes with closed timelike 
curves, and warp drives are forbidden (FSW93tlHaw92j . The WEC is also used in theorems, 
which prove that at a certain stage the formation of singularities becomes inevitable. These 
singularities are related to gravitational collapse |Pen65j as well as to an initial big bang 
[HE73J. Only recently, it was suggested that in inflationary cosmology an initial singu- 
larity may exists even if the energy conditions are violated [B GV03| . Most important for 
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i 



(4.1) 



T av u»u u > 
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us, energy conditions are necessary for the second law of thermodynamics to hold (see the 
discussion below). 

As opposed to classical field theory, it is known that the weak energy condition is violated 
in quantum field theory. In fact, in any Wightman field the renormalized energy density 
may become arbitrarily negative at a spacetime point [EGJ651 IKuo97] . This may happen 
under very general assumptions and for free as well as interacting local quantum fields. 
Simple examples of states with negative energy density are given by certain superpositions 
of the vacuum state with a two-particle state |Pfe98t IFHR02] . Besides this, there are also 
examples for negative energy densities that do not rely on the availability of particles, e.g., 
the Casimir effect or squeezed states of light (see references in |Fla97j ). 

Quantum inequalities were introduced originally in |For78j . where it is argued that un- 
constrained negative energy fluxes F = E/t can violate the second law of thermodynamics. 
However, all evidence tells us that the second law of thermodynamics holds on a macro- 
scopic level. Quantum inequalities provide a mechanism that prevents the microscopical 
violations to survive on a macroscopic level. The argument goes as follows. An absorber 
is a quantum system that has an energy uncertainty AE > £ -1 , where t is some timescale. 
The magnitude of change in its energy due to absorption of negative flux is of the order of 
t\F\. Hence, no macroscopic effects remain if the magnitude and duration of the negative 
energy flux is constrained by \F\ < t -2 . 

To constrain the negative energy that observers can measure it was suggested using the 
integral of the energy density over the worldline of a geodesic observer |Tip78| - a concept 
which led to the averaged weak energy condition (AWEC) dr (T^ u u^u u ) > 0, where u^ 1 
is the observers four-velocity and r is his proper time. This condition allows for the energy 
density to be pointwise negative as long as there is enough positive energy elsewhere on the 
worldline to compensate for this. A related condition is the averaged null energy condition 
(ANEC), which follows from integration along null geodesies [WY91J . Unfortunately, these 
conditions are violated easily for the vacuum state in certain spacetimes. It was then 
discovered [FR95j that introducing a sampling function into the AWEC integral successfully 
constrains the magnitude and duration of negative energy densities. The idea of averaging 
the renormalized stress-energy tensor over regions or curves in spacetime has been very 
fruitful since then and led to a variety of quantum weak energy inequalities of different 
types in different settings - see the reviews [Pfe98t IFew05a| . 

An absolute quantum energy inequality (QEI) for the renormalized stress-energy tensor 
has the general form [FP06J 



where the integral is taken over a region K of spacetime and the sampling tensor belongs 
to the class of second rank contravariant tensor fields, possibly singularly supported on 
curves or surfaces in /C. Furthermore, Q is a positive, real-valued map on these tensors, 




(4.3) 



K 
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and oj belongs to a class of states of the theory - usually the class of Hadamard states if 
curved spacetime is considered. 

Most quantum energy inequalities are obtained by averaging with a sampling function 
fit) along the worldline of an observer, 

J dt fit) (: : (M)>^KW > -Q(f) , (4-4) 

7 

where denotes the four- velocity of the observer. Early quantum inequalities were derived 
for Lorentzian sampling functions ft {t) = t / (rf{t 2 + £§)). To give an explicit example, for 
a free quantized scalar field in Minkowski spacetime it holds 

oo 

I dt f t0 (t) (: T, u :) uV > -3^574 > (4-5) 

—00 

in the frame of an arbitrary inertial observer with proper time coordinate t for any state 
|FR95] . In the meantime, the class of admissible sampling functions has been extended to 
arbitrary smooth positive functions with sufficiently nice decay properties. 

A simple configuration leading to violation of the AWEC and absolute QEIs is an observer 
at rest between uncharged, perfectly conducting plates. Then, the Casimir effect |Ca s48j 
predicts negative vacuum energy for the quantized electromagnetic field measured by such 
an observer. In curved spacetimes a further difficulty is normal ordering of the stress-energy 
tensor. This problem is solved by difference QEIs, which are inequalities of the type 

J dt fit) [(T^(t))) u - (T^mj Y(t)Y(t) > -Q(j>o) , (4.6) 

7 

where ujq is a reference state. Clearly, in the Fock representation of a Hadamard state ujq 
we would, by Wald's fourth axiom (see section [2.2.41) . regain an absolute QEI. Difference 
QEIs have been proved on arbitrary globally hyperbolic spacetimes with general sampling 
functions on the class of Hadamard states. See below for a precise statement. 

Sampling along timelike curves is not the only possibility to obtain energy inequalities. 
One can try sampling in a spacetime region or one may sample over a spatial region alone: 

H(£, E) := J ^(x)/(x) (: T, u :) , (4-7) 

s 

where £ M is a vector field that is orthogonal to S and /(x) has compact support on S. Only 
for a restricted class of theories, namely, conformal fields in two spacetime dimensions, 
spacelike sampled quantum inequalities have been derived [Fla97l IVolOOj . Rather it seems 
that compactly supported weighted averages over spacelike surfaces alone are generally un- 
bounded below for dimensions n > 3 [FHR02J. Naively, this kind of objects could have been 
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hoped to be prototypes of 'local Hamiltonians'. Interestingly, the integrals involved here 
resemble the integrals that were investigated some forty years ago in order to derive sym- 
metry generating global charge operators from local currents (see |V5177t IOrz70t Req76| ) . 
In particular, these investigations suggest, that spacelike smearing alone of local currents, 
does not yield operators with sensible properties. In view of these problems with spacelike 
averaging, we adopt the viewpoint that the investigation of timelike averaged stress-energy 
tensor energy densities is more promising. 

Let us make two more remarks. First, we note that there are quantum field theories 
which do not satisfy quantum inequalities. For example, the non-minimally coupled scalar 
fields violate the energy conditions already on the classical level and, as expected, their 
averaged energy density is unbounded below on the class of Hadamard states. Recently, 
Fewster and Osterbrink derived state dependent quantum inequalities for the fields with 
coupling < £ < \ |FU08j . Second, the issue of quantum energy inequalities for interacting 
fields is not yet resolved, though, recently, some progress in this direction has been reported 
by Fewster and Bostelmann [FBj . 



4.1.1 Stability conditions 

Quantum inequalities seem to be closely related to other stability conditions in quantum 
field theory, as has been pointed out by Fewster and Verch |FV03l IFew05b] . On the micro- 
scopic level the microlocal spectrum condition (theorem 12.6 j) serves as a suitable stability 
condition. Now, theorem 14. II states that quantum inequality exists for all states that satisfy 
the microlocal spectrum condition on arbitrary globally hyperbolic spacetimes. In [FV03j a 
macroscopic stability condition was related to the quantum inequalities. It was shown that 
on static spacetimes the existence of passive states (see below) follows from the existence 
of quantum weak energy inequalities. The circle is closed by noting that in [SVOOj it was 
shown that the two-point distributions of passive states of (vector-valued) quantum fields 
in a stationary-spacetime satisfy the microlocal spectrum condition. Now, viewing quan- 
tum inequalities as mesoscopic stability conditions this relates stability conditions on three 
different scales. 



Passivity 

Let us briefly comment on the notion of passivity. A state on a C*-algebra (21, at) is called 
passive if and only if 

u (u*j(U)j > (4.8) 

for any U G W(2l) fl D(5), where IA Q denotes the unit-component of the group W(2l) of all 
unitary elements of 21 and D(5) is the set of all A 6 21 for which the the generator 

5(A):=hm- t (a t (A)-A) (4.9) 
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exists |PW78[ IBR97] . For example, for bounded operators on the Hilbert space of a 
quantum mechanical system with Hamiltonian H, the generator amounts to a bounded 
symmetric derivation 6(A) = i[H,A] (see the monograph [Sak91] ). The notion of passiv- 
ity was introduced in |PW78j as a precise mathematical formulation of the second law of 
thermodynamics, which says that systems at equilibrium are unable to perform mechanical 
work in cyclic processes. This justifies to call it a macroscopic stability condition. It is 
valid for infinite systems and is closely related to the KMS condition (see chapter [5]) by the 
fact that KMS states and mixtures of KMS are passive. Under certain technical conditions 
ensuring that we deal with pure phases, the inverse statement, namely, that passive states 
are KMS states or ground states, is also true. 

To date, there exists no formulation of passivity for non-stationary spacetimes. A direct 
implementation would require to replace the (strongly-continuous) group of automorphism 
a t by a propagator family of automorphisms a tjS (see the related discussion in section |2.2j) 
and thus to deal with a time-dependent family of derivations 5 t - Such dynamical families 
are not well investigated with respect to (non)-equilibrium states - only two references 
are know to the present author. First, in [BMS02J certain states on Robertson- Walker 
spacetimes, which were obtained by 'transplantation' from de Sitter space, are proved to be 
locally passive in a certain sense. Second, in |Oji86| a general framework for the treatment of 
time-dependent non-equilibrium processes was proposed, which makes use of the propagator 
families and their generators. This (perturbative) scheme is equivalent the standard thermo 
field dynamics if and only if the states under consideration are equilibrium states and the 
dynamics is time-independent. This remark ends our digression on passivity. 



The classical energy-momentum tensor of the minimally coupled massive scalar field is given 



The quantized and point split version of (I4.10p is defined as follows |Few00t IFS07] . For 
a smooth timelike curve 7(i) in Ai, where t is the proper time of the curve, denote by 
U a tubular neighborhood of 7. Choose an orthonormal frame {e^} a= o, 1,2,3 in U so that 
Qfiu _ n a/3 e v e v anc j suc j 1 that the restriction of to 7 equals the four-velocity of the curve, 
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by 



T„ v = V M 0V^ - -g^ ( V ff 0V CT + m 2 2 ) 



(4.10) 



e o I7 



r(t). 
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Suppose that the two-point distribution of too obeys the microlocal spectrum condition. 
Then one can define a distribution by 

1 3 



^0 - 

a=0 



+ -mV* [(1®1)4 2) (^^)] ■ ( 4 - n ) 

Here Lp* : M x. M —> R 2 denotes the pull-back induced by the map <p(t,t') = (j(t), j(tf))- 
That (14.111) is a well defined distribution on R 2 is shown in [FewOOj . Let 21 be the algebra 
of observables of the minimally coupled scalar field of mass m > on a globally hyperbolic 
spacetime {M.,Q^ U ) with dimension n > 2. Let, furthermore, 7 : R — ► M. be a smooth 
timelike curve. Then the following theorem holds. 

Theorem 4.1 ( jFewOO] ) Let u and oj$ be states on 2l(.M, 0^) with globally Hadamard 
two-point distributions and define the normal ordered energy density relative to ujq by 

('■ T ') u := (T) u - (T) Uo . (4.12) 
Then (: T :) is smooth, and the quantum inequality 



OC 



J dt fit) 2 (: T :> u (t,t)>- | ^ [(/ ® /) <T> J (-A, A) (4.13) 

7 

ZioZds for all real-valued f G @(R), and the right hand side is convergent for all such f . 

4.2.1 Stress-energy in Robertson-Walker spacetimes 

As a necessary prerequisite for the construction the almost equilibrium states we need the 
explicit form of the left hand side of (I4.13P in Robertson- Walker spacetimes; This will be 
part of our free energy functional to be defined in chapter [51 We state the result as a lemma. 



Lemma 4.2 Let T^ u be the stress-energy tensor of real, linear, scalar field in Robertson- 
Walker spacetimes and let uj^ujq two homogeneous, isotropic states. Furthermore, let f e 
3>(M) and 7 be the worldline of an isotropic observer. Then, 

J dt f(t) 2 ((T) u (t, t) - (T)^ (t, t)) = J d^kf J dt fit) 2 ( Pk (t) - p^it)) , (4.14) 

where 

P k (t) :=h{S) (\T k (t)\ 2 + ut\T k (t)\ 2 ) +^{b 2 (S) (T k (t) 2 + ^T k (t) 2 )} , (4.15) 
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and the coefficients 61,62 are defined in (I3.61al) . The energy density p^ (t) is defined anal- 
ogously. The measures dfi(k) £ differ by constants from the measures dk defined in \3.21\ 

Hk dk dk 

Mk) ° ■= ' Mk)+ := W + ' ^ " W • (416) 

Proof. The stress-energy tensor of the Klein-Gordon field (12. 2p is given by equation (14.101) . 
In order to calculate the left hand side of (14.141) we have to consider the point-split expression 



J dt f(tf [(T m )Jt,t) - (T x ) a Jt,t)] 

= /* /w *ita(iv v;-i£ 



(4.17) 



x=-y(t) 
x'= 7 (t') 



where 

F(x,x') :=cu (2 \x,x') -Jq\x,x') (4.18) 

and the two-point distribution cu^(x, x') is given in lemma 1X21 We need the limits x' — > x 
of the involved derivatives of the functional F(x,x'). For the restricted case of pure states, 
this calculation has been performed in [Olb0 7b|. See also |PF74] for similar results. In 
section 15.51 we will see that the limiting case of pure states amounts to the values 62 = 
and 61 = |. 

The difference F( smooth function if oj and lu are Hadamard states. In that 

case, the coincidence limit x — > x' is well defined. The derivatives involved in the energy 
momentum tensor and which have to be calculated for each of the three cases e = 0, 1, — 1 
are 

3 

lim F(x,x') , lim V V F{x, x') , lim ^ V M V' t/ F(x, x') . (4.19) 

With the convention that latin indices are summed from 1 to 3 we can simplify the spatial 
derivative as 

3 1 

rv,v: = f>tfv f v; = . (4.20) 

a 

H,v=X 



Besides that, we simplify the notation by introducing ellipsis for the second summand, 
which always looks like the first except that it belongs to the state ujq. 
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[e = 0]: This is, as expected, the simplest case. The modulus |Yk(x)| 2 = -^p- of the 
spatial solutions is independent of x. Thus, we can calculate directly 



Ux°)\ 2 { h + ^ ) + T k (x ) 2 ■ b 2 



lim F(x,x') = [ dk |y k (x) 

x'— *x J 

+ T^°Y-h + \Tv(x°)\ 2 (b 1 

dk [2\T k (x°)\ 2 ■ h + {T k (x ) 2 ■ b 2 }} 



(2tt)' 



2 1 d/x(k) [|T k (x°)| 2 • 6! + U (T k (x ) 2 • 6 2 }] - 



(4.21) 



(4.22) 
(4.23) 



where we defined the measure d/x(k) 
time derivatives 



dk 



Similarly we obtain for the term involving 



lim V V F(x,x') = 2 / dfi(k) \T k (x°)\ 2 ■b l + ^\ T k (x ) 2 ■ b 2 



(4.24) 



The third object involves spatial derivatives. Is is calculated with the aid of f)4.20p as 



lim V q^V^V u F(x,x' 

x'^x ' * 

1 



— — lim s'NiV.F^x') 
a(x u ) / x'^x J 



ax 



0\2 



dk s^V l F k (x)V;r k (x') [2|T k (x )| 2 • h + 2K {T k (x ) 2 ■ b 2 }] 



(27F / dk aW [2|Tk( "° )|2 ' hl + m {TkiX ° )2 ' h}] ~ 
2 / [\Ux°)\ 2 • 6i +3? (T k (x ) 2 • b 2 }] - . . . . 



(4.25) 

(4.26) 
(4.27) 
(4.28) 



Here we used that s ij 'V i F k (x)V j F k (x) = ^ • E(k), which follows directly from (13321 . 



[e — +1]: In this case we have to work a little more. The states under consideration are as- 
sumed to be homogeneous. Thus, the quantities lmv^ F(x, x') and \\m x i^ x VoVg-F(a;, x') 
cannot depend on the spatial coordinates x and x'. By the compactness of S + , we can 
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carry out the limiting procedure, integrate over S + using the measure dx and divide by the 
volume Ve+. This gives 



lim F(x, x') 



1 



Vy. 



dx lim Fix, x') 



dk / dx \Y^)\ 2 \\T k {x°)\%+T k (x»)% 



e+ 



+ T k (x°)% + \T k (x°)\%] - 

= / dk5(k,k) [2|T k (x )| 2 6 1 + 2K{T k (x ) 2 6 2 }] 
- 2 / rf/i(k) [\T k {aP)\% + to {T k (x°)%}} 



(4.29) 
(4.30) 

(4.31) 
(4.32) 



where we used (13.291) and defined the measure dfi(k) 



dk 



Again, after taking the time 



derivatives, we can perform the same calculation as before: 

Jim V V F(x,x') = 2 J dfiQc) [\f k (x°)\ 2 b 1 + n\f k (x°) 2 b 2 j 
The spatial part gives the preliminary expression 

lim V tf v V^' v F{x,x') = J- / dx— j^- lim 5 ^V l V'F(x,x / ) 

s'^x * v v Ve+ / a 2 (x°) x'^x 



(4.33) 



1 1 

a(a; ) 2 



dx / dk lim s ij Vilk(x) V'F k (x') 



(4.34) 
(4.35) 



E+ 



X 



1 



T k (x°)T k (x'°) [&! + -)+ T k (x°)T k (x'> 2 



+ T k (x°)T k (x'°)b 2 + T k (x°)T k (x'°) [ h 



1 1 



a(x ) 2 



dk / dx s 11 V i Y k (x)V j y k (x) 



(4.36) 



x 



iTk^ )! 2 ^ +T k (x°) 2 6 2 + T k (x°) 2 6 2 + llUx )!^! 



Since S + is compact without boundary, we can perform a partial integration |Tay96| using 



dx s v V i Y k (x) V,r k (x) = - J dx AF k (x) F k (x) . 

E+ E+ 



(4.37) 
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Furthermore, since A s Yk(x 

1 ' dk J dx |r k (x 



1 

2 



-£'(k)Y k (x), it follows that 
|2 E(k) 



a(x ) 2 



2\T k (x°)\ 2 b 1 + 2^{T k (x°) 2 b 2 } 



dk <5(k, k) 
£(k) 



a(:r ) 2 



goo 

&(x ) 2 
\T k (x 



2\T k (x°)\ 2 b 1 + 2^{T k (x°) 2 b 2 } 



0m2 



6 1 + 3?{T k (x°) 2 6 2 } 



(4.38) 

(4.39) 
(4.40) 



[e — — 1]: In this case one uses techniques from the harmonic analysis of hyperbolic 
spaces with negative constant curvature, also called Lobachevskian spaces |GGV66j . In 
order to find a direct integral representation of the isometry group G~ , one embeds X - as 
a three-dimensional hyperboloid with coordinates x = (xi, x 2 , £3) into the four- dimensional 
Minkowski spacetime (M, n^) by means of the map 1 : R 3 — > M, x 1— ► + x 2 , x). One can 
calculate the metric on £~ as the pullback g M „ = L*n^ v , which gives 



M x ) 



ff J (x) = <L- + XjX 



(4.41) 



Now, define a normalized momentum by £ := k/k and write for each k the Fourier transform 
of a function h e £^(E~) as 

h k (£) := y dx F k€ (x)/i(x) , (4.42) 

where lk(x) are eigenf unctions of the Laplace operator given by (13.341) . This Fourier trans- 
form is a map with values in L 2 (S 2 ,dfl), where S 2 is the two-sphere embedded in R 3 and 
dQ denotes the induced measure on S 2 . The Lorentz transformations g : f(x) 1— > f(xg) on 
the functions f(x) correspond to operators 



U k (g)h k = [U(g)h 



(4.43) 



Let g x be the Lorentz transformation that maps (vT+x, x) to (1,0,0,0). Then by the 
unitarity of U k (g) 



\YJx) 



(2tt) 3 
1 

(X) 3 
1 

1 



dm {x-g-^y 

dfi(£) (g x x • 

47T 



(4.44) 
(4.45) 
(4.46) 
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Again, as for e — 0, the modulus of lk(x) is a constant. Using this, we obtain 
lim F(x,x') = [ dk |F k (x)| 2 



„0\2r 





-T k 








-1) 







An 



dk [2|T k (a; )| 2 + 2K{T k (a; ) 2 &2}] 



2 / dfiQs) [\T k (x°)\ 2 + 3? {T k (x°) 2 b 2 }] 



where d/x(k) = Similarly, one finds 

lim V V' F(x, x') = 2 / d/x(k) |T k (x°)| 2 + 3? <| T k (a; u ) 2 6 2 



Finally, since s^V i Y k (x)V'y k (x / ) = (1 + k 2 )|F k (x)| = £(k)|F k (x)| we obtain by 



YJx) 



1 / o x-| 

x — 



(2tt)5 V k 
(see (13.341) ) for the spatial derivatives the expression 

3 



-l+ik 



lim y grV^Ffa x' 

x'— *x z — 4 
1 



a 2 (x°) z'^z 



lim s ij ' Vi VjF(x, x') 



£i[ X 



0\2 



dk ^v l y k (x)v;.y k (x / ) [2|r k (x°)| 2 + 23? {T k (x°) 2 b 2 }} 



(27F / dk aW [ 2|Tk(x ° )|261 + 2K i T ^°) 2&2 }] - 

2 /rfMk)4nKO r ^ )l 2fo i + ^{ r k(x ) 2 &2}]-... • 

J a [x ) 



(4.48) 



(4.49) 
(4.50) 

(4.51) 



(4.52) 



(4.53) 

(4.54) 

(4.55) 

(4.56) 
□ 



Putting all this together gives the desired result. 
A better parametrization 

For the actual minimization, it will be convenient to parametrize the averaged stress-energy 
by the components of the two-point matrix S. For this, note that 



dt f{tfp k {t) = hiSMT^f) + K{& 2 (£)c 2 (T k ,/)} 



(4.57) 
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where we have defined 

Ci(T k ,f) := J dt f{tf (\f k (t)\ 2 + u 2 k \T k (t)\ 2 ) j (4 . 58a) 
c 2 (T k , /) := J dt f(tf (f k (t) 2 + u 2 k T k (t) 2 ) . (4.58b) 

Our definition of the quantities ci(T k , f) and c 2 (T k , /) differ from the one in [Olb07a] by a 
factor of 2. Now, inserting b\ and b 2 we can write 

J dt f(t) 2 p k (t) = S 00 di(T k , /) + a 6 5 n d 2 (T k , /) + a 3 S* d 3 (T k , f) , (4.59) 

where 

d 1 {T k ,f) := \T k (t )\ 2 Cl {T k J) -&{Tl(t )c 2 (T k J)} , (4.60a) 

d 2 (T k J) := |T k (t )| 2 Cl (T k ,/) -^{^(t )c 2 (r k ,/)} , (4.60b) 

rf 3 (T k , /) := (T k (t )f k (t ) + T k (t )t k (t )) ci(T k , /) - 3? {2T k (t )t k (t )c 2 (T k , /)} . 

(4.60c) 



5 Almost equilibrium states 



In this section we tackle the construction the almost equilibrium states. First, we explain the 
basic ideas behind the procedure and the definition of the free energy functional associated 
to the states of interest. Then, we calculate the entropy, which is the part of the free energy, 
that cannot be inferred from the quantum energy inequalities discussed in the last chapter. 
This is followed by the actual minimization and the definition of the almost equilibrium 
states. In the last section we prove that the almost equilibrium states are indeed Hadamard 
states. 



5.1 KMS states 

Consider a finite quantum system that is described by an algebra of observables 21. A 
general state uj on 21 is described by a density matrix p and the expectation value of an 
observable A e 21 is given by 

u(A) = Tr(pA) (5.1) 

Now, a state u on 21 is a thermal equilibrium state at inverse temperature (3 G R if it is a 
Gibbs state. A Gibbs state is given by a density matrix of the form 

e~ 0H . ■ 

PpH : ~ Tr(e-^) ' ^'^ 

where H = H* is a self-adjoint operator (the free Hamiltonian) and e _/3H is a trace class 
operator, which means that Tr(e _/3H ) < oo. The trace class property of e~@ H is guaranteed 
for finite volume systems by the properties of H in that case [Haa96] . A simple calculation 
using the cyclic invariance of the trace, Tr(ABC) = Tr(CAB), shows that Gibbs states 
formally satisfy the relation 

u(AB) = ou(Ba ip (A)) , (5.3) 
for all A, B 6 21, where the automorphism 

a t (A) := e itH Ae-' ltH (5.4) 

describes the free time evolution of an observable A e 21. Equation f)5.3p gives the combina- 
torics inherent to thermal equilibrium states of finite as well as infinite quantum systems. 
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Now, consider a state u p that is described by a density matrix p and define the entropy 
of Up by the von Neumann entropy functional 

S(u p ) :=-Tr(plnp) , (5.5) 

where we understand that xlnx = at x = 0. If one defines the free energy of the state uj p 
by 

T{u p ) := oj p (H) - , (5.6) 

then the Gibbs state is the unique state that minimizes the free energy T{oj p ) . Equivalently, 
the Gibbs state maximizes the entropy at fixed energy |BR97, Weh78j . 

A necessary condition for Gibbs states to be well defined is that H has a purely discrete 
spectrum bounded from below (We consider only positive temperature states: (3 > 0). 
For infinite systems, Gibbs states are ill defined because Tr(e - ^) = oo. Nonetheless, a 
generalized definition of thermal equilibrium states exists for infinite systems, namely, the 
notion of KMS states. The above notion of a Gibbs state implies a characterization in terms 
of analytic functions F Ai b f° r an Y pah A, B £ 21, that satisfy certain boundary conditions, 
provided that the Hamiltonian H is bounded below and e _/3H//2 is of trace class (see, e.g., 
|Emc72t IWeh78j ). That this characterization survives the thermodynamical limit and thus 
can be used for the definition of equilibrium states in the general case has been shown 
in [HHW67J. Since then, many applications have shown that KMS states indeed describe 
thermal equilibrium. 

We come to the definition of KMS states. Consider a C*-dynamical system (21, at). 

Definition 5.1 (KMS state) Let (3 > denote the inverse temperature. A state up is 
called a /3-KMS state if for all A, B £ 21 there is a function F a ,b{^) '■ Cp — » C which is 
analytic in the strip Cp := {z £ C : < ^sz < (3}, bounded and continuous on its closure 
Cp, and satisfies the boundary conditions 

F A>B {t) = up{Aa t {B)), F A>B (t + i(3) = u p {a t {B)A) . (5.7) 

Of course, on a non-stationary spacetime there exists no global C*-dynamical system. Thus, 
no global KMS state can exist in that spacetimes. We remark that a relativistic KMS 
condition has been defined in H BO II . 

5.2 Almost equilibrium states 

As noted in chapter HI it is well known that in non-stationary spacetimes there is no global 
generator of time-translations. The standard Hamiltonian - the integral of the energy 
density over the Cauchy surface - is not conserved. Even worse, the energy density is 
not bounded-below at a spacetime point, and spatially averaged energy densities are ill 
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defined for dimensions n > 2. An energy quantity which is unbounded below, poses serious 
problems regarding the existence of equilibrium states. However, as we have also seen, there 
exists a quantum inequality that gives lower bounds on the difference of energy densities 
of Hadamard states averaged along the worldline of an observer in any globally hyperbolic 
spacetime. Thus, one can expect to find ground states with respect to this kind of energy 
Indeed, in |Olb07a] states with this property so-called states of low energy, were constructed 
(see section 1531) . However, states of low energy are pure states; consequently, they do not 
describe systems in thermal equilibrium. They are merely ground states of a class of almost 
equilibrium states, as we will show in this chapter. 

Our guiding idea for the construction of the almost equilibrium states is to define a sensible 
free energy functional J-{u) on the class of homogeneous, isotropic, quasi-free states that is 
based on the quantum inequality stated in theorem !4.13l We write this functional exclusively 
in terms of the two-point matrix (13.491) . Subsequently, we minimize this functional with 
respect to the components of (13.491) . 

Suppose that we can define a sensible free energy functional J-"(u)), where we replace the 
energy expectation value uj p (H) present in (15. 6p by the worldline averaged expectation value 
of the (renormalized) stress-energy tensor, 



for some real- valued / G @(R). Additionally, we assume that 7 is the worldline of an 
isotropic observer. To remind the reader, an isotropic observer moves along a timelike 
geodesic with tangent vector orthogonal to the homogeneous spacelike surfaces. For such 
an observer, the spacetime looks spatially isotropic at any instant of time. Note also that 
T{oj) is actually a free energy density rather than a free energy and, at this stage, it is 
unclear how to define the entropy S(u). 

The minimization of F{uj) means the following. We consider the class of homogeneous, 
isotropic, quasi-free states. All such states are completely determined by their two-point 
distributions (13.601) . Within this class we look for a distinguished state ujpj such that for 
any other state uj in the same class the free energy is larger, i.e., 



In the rest of this chapter we will show that such a state up j exists and is uniquely deter- 
mined by its two-point distribution. Furthermore, we show that the two-point distribution 
of the state uipj satisfies the Hadamard condition. This latter step is crucial for the whole 
construction by the following argument. 

The difference of the energy densities involved in the left hand side of (15. 9p is 




(5.8) 



7 




(5.9) 



ou(H f ) - u> M (H f )) = Jdt f{tf (<T> W (M) - (T) UfiJ (M)) 



(5.10) 



7 
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Provided that the states u> and Upj satisfy the Hadamard condition and dim.M > 2, the 
integral on the right hand side of (I5.10p is well defined and bounded from below as we know 
from theorem (I4.13p . At this stage, it IS, clS db matter of fact, not known if the state that 
eventually minimizes this functional is a Hadamard state. So, the initial ansatz is justified 
afterwards by proving that the resulting state Upj is indeed a Hadamard state. 

The integral on the right hand side of (15.101) has the form given in lemma 14.21 It is 
a difference of integrals of the energy density p k for each mode k. Note that this mode 
decomposition of the energy is due to the symmetry of the spacelike surfaces in Robertson- 
Walker spacetimes and familiar from standard quantum field theory, which, after all, can 
be viewed as a theory of infinitely many, independent, harmonic oscillators. We infer from 
this observation that, when looking for the minimum of the averaged energy density, we 
can confine ourselves to a single mode k. Put differently, we can minimize the integrand 
instead of the integral to obtain the overall minimum. In order to utilize this principle for 
the free energy, we write 

:= J d/i(k)^ k , (5.11) 
where we have defined the free energy of a single mode by 

■= J dt f(t) 2 p k (t) - i<S k (u;) . (5.12) 

7 

We have yet to define the entropy 5 k . Since each mode k can be considered as an indepen- 
dent quantum system with one degree of freedom, it is natural to ascribe an entropy given 
by the von Neumann entropy functional (15.51) to it. Then, 5 k is completely characterized 
by the density matrix one ascribes to the single-mode system. Since we want our states to 
exhibit thermal behavior, we take our modes to be Gibbs equilibrium states. This means, 
we consider density matrices 

where we assume that K is a positive definite, quadratic form that is not diagonalized from 
the outset. The latter point is important, since, such a form is determined by three real 
parameters, just as the two-point distribution (13.491) of a quasi-free state. 

The further course of action will be as follows. In lemma 15.21 we state the generating 
functional of KMS states with respect to the evolution generated by K. Then (by equation 
(I5.54p ) we obtain a one-to-one correspondence between the generator K and the two-point 
matrix S. This relation is used to express the free energy jF k in terms of Sii, 5*22, and [S], 
which then allows to minimize with respect to these variables. The result of the minimization 
will be a uniquely determined state of inverse temperature (3 associated to the sampling 
function fit). 
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5.3 The entropy 

For the entropy of a mode k, we use the von Neumann entropy associated to the Gibbs 
state with density matrix p^x- 

S\,(uj pf)K ) = -Tr(p m lnp m ) , (5.14) 

where we have defined a positive definite, real quadratic form for the position and momen- 
tum operators (q,p) by 

K(q,p) := K 00 q 2 + K 01 (qp + pq) + K 11 p 2 = (£™ • Q • (5-15) 

In order to calculate the entropy we use canonical diagonalization of K. This is possible 
for any positive definite quadratic form on the phase space. 



Canonical diagonalization of K 



The canonical diagonalization of a quadratic Hamiltonian is a well investigated technique 
(for general results see, e.g., |MQ71[ IBG79] ). We consider the simple case of linear trans- 
formation matrices M for the coordinates q,p that are canonical, which means 



M J JM = J 
The eigenvalues of the matrix 



J 



1 
-1 



JK 



K i 
-K 00 



-K m 



(5.16) 



(5.17) 



are ±i^/[K]. It turns out that in this model Vt := 2y[iT] plays the role of the oscillator 
frequency. The eigenvectors v = ni] + im 2 and v = ni! — im 2 are 



-K 01 







-K m 



(5.18) 



which defines the vectors mi, m 2 needed for the transformation matrix M = (mi, m 2 ): 



M 



y 



K n 
-K 01 







The normalization constant y is fixed by the condition mJJnij = Jy, which gives 

1 

y = z 



K 



ii 



(5.19) 



(5.20) 
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so that 

M = ± — ! ( K £ Lr) ■ (-21) 

This M diagonalizes the quadratic form K by 



M T KM = {^f JL- } ) . (5.22) 



The new coordinates are given by 

1 



Q=^^r^yr^v^i^, (5. 23) 



" ■ 1 f °Tf P> ) (5.24) 



If we define the frequency : = 2 y^if] then the diagonalized matrix reads § ( q n ) ■ How- 
ever there is still some freedom left in the choice of the coordinates. We can perform a 
simultaneous transformation q' := yQ and p' := y~ l P in order to gain a more suitable form 
of the matrix. Choosing y = y/Q leads to 

K = J °) • (5.25) 



2 V° 

We conclude that choosing appropriate coordinates (q',p') is is possible to diagonalize 
the 'Hamiltonian' matrix K and write it in the standard harmonic oscillator form 

K( q ',p') = ~(tfq' 2 +p' 2 ) . (5.26) 

Therefore, for all basis-independent quantities, we can use standard results from the the 
theory of the one dimensional, harmonic oscillator with frequency Q. For example the 
eigenvalues of the system are 

K n =(n+~^2y/[K], (5.27) 
and the partition function is (see, e.g., [Rei87j) 



P -P\/IK] i 

Z := Tr(e^ K ) = — = = ? r- . (5.28) 

1 _ e -Wy/lK\ 2smh((3 x /\K]) 



From the partition function Z one obtains the entropy as 



S k (u> P0K ) = - In (l - e - 2 ^) + 2/3 y/\K\- 



(5.29) 
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5.3.1 The generator of KMS states 

Now we know the entropy S^{uOp„ K ) of the modes in terms of the Hamiltonian matrix K. 
Note, that the expression (15.291) involves only the the determinant of K. Next, as we want 
to characterize states by their two-point distribution, we need the correspondence between 
the matrices K and S. For this, we characterize the KMS states on a Weyl algebra with 
respect to the time-evolution generated by K. The result, which is stated in the following 
lemma, generalizes a result of [NT93] . 

Lemma 5.2 Let 21 be the Weyl algebra generated by the exponentials of the position and 
momentum operators and let the time evolution be generated by a positive Hermitian form 
K, with [K] := detiT ^ on the phase space elements z e R 2 . Then, the generator of a 
KMS state associated to the inverse temperature (3 > on 21 is given by 



u(W(z)) = e~^ [K] zTK ~ lz cot HV WP) . (5.30) 

Proof. To prove this, we consider the exponentiated one-dimensional Heisenberg *-algebra 
generated by q and p, which are subject to the commutation relations 

[q,p]=i, (5.31) 
[q,q] = 0=[p,p] , (5.32) 

and invariant under involution, q* = q, p* = p. The Weyl operators are given by 

W(z) := e - i{ziq+Z2p) , (5.33) 

where z = (zi, z 2 ) G R 2 . We define a symplectic form a by 

cr(z, z') = z\z' 2 — z 2 z 1 . (5.34) 

The Weyl operators W(z), which generate the algebra 2U, satisfy the Weyl relations 

W(z)* = W(-z) , (5.35) 
W{z)W{z') = e'^ (z ' z,) W{z + z') . (5.36) 

A state u is, as usual, defined as a positive, linear, normalized functional on 2U. By 
the GNS theorem we have a representation 7r w and a representation space 7i w with vectors 
|z) = ty(z)|0) and |0) being the cyclic vector of TC U . A state u on 2U is determined by its 
action on the Weyl operators 

f(z):=u(W(z)) = (0,W(z)0) . (5.37) 
The scalar product of two arbitrary vectors is given by 

(z', z) = (W(z')0, W{z)0) = e^ {z '- z) /(z - z') . (5.38) 
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The time automorphisms at over 2B correspond to symplectic transformation matrices B t 
on the phase space. The KMS condition (15.71) together with the Weyl relations yields 

e -i CT (z,z')j( z + z /j = e -H*'< B ^)f(B if3 z + z') . (5.39) 
If z' = — z this simplifies to (Note that /(0) = 1): 

1 = e^ z ' B ^ z) f((B if3 - l)z) (5.40) 

f((B if} - l)z) = _ ^ 

Since for invertible (5^ — 1) we have 
a transformation (5^ — l)z — > z gives 



-77 cr 



/(z) =e 2 V^- 1 ; . (5.43) 

Next, we tackle the relation between the matrices .fT and S. The time evolution generated 
by K is given by 

a t (W(z)) = e -^t e i(z iq +z 2P ) e iKt = e i(zi(t) q +z 2 (t) P ) _ ^ 5>44 ) 

The time evolution of z(t) is given by the Heisenberg equation z = i[K, z]. The commutators 
[K, q) and [i^p] can be calculated by flSTBTD and the identity [A, BC] = [A, B]C + B[A, C). 
The Heisenberg equation gives the system of differential equations 



zi\ = (ZK m -2K 00 \ (z x 
z 2 \2K n -2K 01 ' U 



(5.45) 



We introduce the frequency 

Q := 2^\K] , (5.46) 

where [K] := K Q0 Ku — Kq X is the determinant of the matrix K. Furthermore, we set 
K i Kqo K n 

K o 1 := 7 =, ^=7^ K " := 7Pr ( ' 

Then, the time evolution can be written as 

z(t) = B t z(0) , (5.48) 

where 

B _ /cos(Qt) + «oi sin(fit) -K 00 sin(fit) \ , , 

4 ' I Kn sin(fit) cos(fit) — «oi sin(fit) J 
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The matrix B t generalizes the usual harmonic oscillator evolution matrix. To implement 
the KMS condition we complexify time, t — > i/3, and since sin(i/3) = isinh(/3) and cos(i/3) = 
cosh(/3), we obtain 

/cosh(fi/3) + i kqi sinh(fi/5) -i k 00 sinh(fty3) \ , , 

i/3 ~ I i k u sinh(fi/5) cosh(^) - i «oi sinh(^) J ' 1 j 



The matrix (.B^ — 1) is invertible for [K] , /3 ^ (note that coth (|) 

i _ 1 A - i Koi coth (^) i K 00 coth (^) 



sinh(z) - 
cosh(x') — 1 < 



(A/3-1)- --- ^ _ i/?iiC0th (M) l +i K 01 coth(f)) ' ^ 

and thus we find 

1 ,zU44=H T f*' "M C !') coth Mxlfl ) ■ (5-52) 



Bip-1 ' / 2 v /[kJ V^^oi ^oo / \ z 2 

So we have found the generator 

/(z) = zTA '" z coth (v») . (5.53) 

This proves lemma f5T2l □ 
Defining the quadratic form S(z,z) := |(z T Sz) and recalling that the generator of a 
quasi-free state is given by /(z) = e~5 5 ( z > z ) we infer that 

(I: fy^m^(^)(-K m (5-4) 

Inverting (I5.54p gives 

(t S^^p^tt "3- (5 - 55) 

The determinants of K and S* are related by 

[K] = ^Arcoth(2v^S]) , vT^] = ^coth(/3 y/[K}) . (5.56) 



Thus, we can write the entropy per mode as 

-2Arcoth(2 v /[5]) 



S*{u PfjK ) = -In (l - e - 2Arcoth ( 2 V^)) + 2 Arcoth(2 — 



e ~2Arcoth(2 x /[S]) 



(5.57) 



= 2 v ^]Arcoth(2 v ^]) + iln(i^L^^) , (5.5S 
where, for the second line, we used that Arcoth(y) = | In f^zfl an d thus e - 2Arcoth (%) : 
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5.4 Minimization of the free energy 

Using the energy density (14.591) and the entropy formula (15.581) we can write the free energy 
functional JF k exclusively in terms of the two-point matrix (I3.49j) . 

Tv ■= S 00 d 1 (T k , /) + a 6 S llC / 2 (T k , /) ± aVSooSn _ [S]d 3 (T k , f) (5.59) 
1 



P 



2^/[S] Arcoth(2v^S]) + i In ' 4 ^ ' 



2 V 4 



where we used [S] := S 00 Su — (Sj^) 2 to eliminate in favour of the determinant [S]. 
The remaining task is a simple exercise, namely, the minimization of this expression. A 
necessary condition for the minimization of JF k is that we have a critical point, i.e., the 
gradient of the function vanishes: 

V ( 5 00 ,5 11 ,[S])- :F k = . (5.60) 

This amounts to the equations 

diiTkJ)+ — 2 — ^s 00 s n -[sr > (5 - 61) 

rf 2 ( Tk , /) + d3{Tk \{ ] = o , (5.62) 

1 ^ 2a(t ) 3 y/SnSn - [S] 1 ' 

1 Arcoth(2 y^I) ± & M 3d f*> D. 1 = Q , (5.63) 



Py/\S\ 2 a/ S 00 S n - [oj 

where the coefficients cfj(T, /) are given by (14.60)) . Multiplying (15.611) by (15.621) we obtain 

\d 



y/SwSn - [S] = 131 J[S]. (5.64) 
A/4aiGt 2 — «3 



Using this, and assuming /3 > 0, the critical point is obtained from (15.63)) as 



VW) = \ COth ^(tfyfedrfk - d^j 



(5.65) 



where we omitted the combination of signs from (15.63)) and d 3 /\d 3 \ giving no solution for 
(3 > 0. So, given a sampling function f(t) there is for every (3 > a unique solution to 
equation ( 15.651) and thus a unique minimum of the free energy functional T{uJ Pf3K ) given in 
(I5.59p . We state this result theorem. 

Theorem 5.3 Consider the Weyl algebra 21 of the minimally coupled Klein- Gordon field in 
a Roberts on- Walker spacetime. Define the free energy of a homogeneous, isotropic, quasi- 
free state u on 21 by 

:= J d/i{k)^ k , (5.66) 
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where 



I 



dt f(t) 2 Pk (t) - • 



(5.67) 



7 



and p k is the energy density of a mode. Averaging takes place along a timelike curve 7 with 
a sampling function f G @(R), and each mode is assigned the von Neumann entropy Sk(u). 
Then, for every (3 > and every function f , there is a unique state upj that minimizes 
F{uj). This state is determined by 



supplemented by equations (I5.6ip . and (I5.62p . where T k is an arbitrary solution to (I3.35P 
satisfying (I3.37p . We refer to LOpj as an almost equilibrium state of inverse temperature [3 
associated to f . 

The almost equilibrium states are parametrized by arbitrary solutions T^(t) of the time- 
dependent part of the Klein-Gordon equation. As we will see in the next section, they 
resemble equilibrium states in Minkowski space if one uses the states of low energy for the 
parametrization. 



In this section, we prove that the two-point distribution of the almost equilibrium states 
satisfies the Hadamard condition. This is of twofold importance. First, it is a necessary 
condition to make our ansatz meaningful, as otherwise the energy difference involved in the 
free energy (15.101) is not bounded below and our formal result remains formal. Second, it 
proves that the almost equilibrium states can be considered as physical states in the sense 
we described in chapter [2J 

We are in the fortunate situation that we can base our proof on previous related work 
|Jun96l [JS021 IUlb07al IOlb07bj . It was shown in [Olb07a] that the states of low energy are 
Hadamard states. The proof used the fact that for large momenta k and large iteration 
order n the difference of the two-point distributions of the states of low energy and the 
adiabatic vacuum states converges to zero. Thus, the two-point distributions differ only by 
a smooth function. Since adiabatic vacuum states of infinite order are Hadamard states 
|Jun96l IJun02j . this proves the Hadamard property of the states of low energy. Along the 
same line of reasoning, we show that the two-point distributions of the almost equilibrium 
states and of the states of low energy also differ only by a smooth function. The calculations 
are simplified a great deal, due to the fact that the states of low energy are the natural 
ground states associated to the almost equilibrium states. 




(5.68) 
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Parametrization by states of low energy 

Let us first introduce the states of low energy |Olb07aj . Consider the Weyl algebra 21 of 
the free Klein-Gordon field over a Robertson- Walker spacetime (M,^) and fit) G 2t(R). 
In the set of homogeneous, isotropic, pure quasi-free states on 21, there is a state oJooj f° r 
which the averaged energy density 

u{H f ) = Jdt f[tf (: T :) u (t) (5.69) 

7 

is minimal. The averaging is understood along the path of an isotropic observer. This state 
is given by the two-point distribution 

<"S/fo*') = / dk y k (x)F k (x / )I k (t)L k (t / ) (5.70) 

with 

L k (t) = \T k (t)+fiT k (t), (5.71) 

where T k is an arbitrary solution of the differential equation (13.351) fulfilling the continuity 
condition ( 13.371) and A, // are given by the equations 



A=exp(!a) VvraF + 5' " = vW^~5' (572) 

The coefficients c\ = ci(T k , f) and c 2 = C2(T k , f), up to a factor of 2, are defined in equation 
(14.58)) . Setting T k = L k , it follows that a state given by L k is a state of low energy if and 
only if it satisfies c 2 (L k , /) = 0. 

For the proof of the Hadamard property, and presumably for most other purposes, the 
expression (15.68)) defining almost equilibrium states is simplified considerably by using the 
states of low energy for the parametrization. So, we plug in a state of low energy L k (t) into 
(15.65)) . Note that the squared modulus of (13.371) . is given by the useful relation 

/— • — \ 2 — • _ 1 

4|T| 2 |T| 2 — [TT + TT] =\TT-TT\ 2 = — j (5.73) 



which turns up several times in the calculations. Plugging in the states of low energy into 
the expression 4d\d2 — d^ gives, due to C2(L k (t), f) = in equation ( 14.601) . 

Ad,d 2 -dj = A\L k (t )\ 2 \L k (t )\ 2 cl(L k ,f) - (l k (t )L k (t ) +L k (toA(to)) 2 c?(^k,/) 

(5.74) 

' cl(L k ,f). (5.75) 



a(£ 
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Using this, equation (15.651) becomes 



[£] = JcotbY ^ Cl (L k J) ) . (5.76) 



2 V 2 

Solving for Sqq and Su at the minimum of the free energy gives 



S u = 7 3M ^/[S] , (5.77) 

a?V Adid 2 — dl 



2a s \d : 



Soo= ,;;r 2 ' = V[S], (5.78) 

where d\ and d 2 have the same sign. Noting that di and d 2 are always positive for states of 
low energy, this reduces to 



1 S 11 = 2|L k (t )| 2 v / P] , (5.79) 
5 00 = 2a(t ) 6 |4(to)| 2 v / M (5.80) 



for the almost equilibrium states. 

Now, we calculate the coefficients 61 and b 2 . We can use the expressions for Soo and Su 
to find an expression of S^ from [S] = Sq Sh — (S 1 ^) 2 : 



5« = ± A /4a(t ) 6 |^k(to)| 2 |^k(to)| 2 - 1 • VIS] . (5.81) 

On the other hand we know that with the normalization (13.371) we have equation (I5.73P 
giving 

a(t ) 3 (L k (t )L k (t ) + L k (i )t k (t )) = ±\/ Mto) 6 \L k (t )\ 2 \L k (t )\ 2 - I • (5.82) 

Each of the last two equations has two possible signs. If the signs adjust such that their 
product is negative we can write 



h = S 0Q \L k (t )\ 2 + a,(t ) 6 S 11 \L k (t )\ 2 + (L k (t )L k (t ) + L k (t )L k (t )j a(t ) 3 ^ (5.83) 
= 4a(t ) 6 |^ k (to)| 2 |4(to)| 2 v / M- (4a(t ) 6 |L k (t )| 2 |L k (t )| 2 -l) VW] (5.84) 



[S] . (5.85) 

Assuming the same adjustment of signs as before, we have an alternative expression for S^, 
namely 



Sg = -a(to) 3 (L k (t )4(t ) + L k {t )L k {t )) y/\S] , (5.86) 
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which inserted into 62 gives 



• 2 



-S 00 L k (t ) - &(to) b S n L k (to) - 2L k (t )L k (t )a(t yS { 



3 0.% 
01 



= 2 a (t ) 6 v^S] \L k (t )\ 2 Ll(t ) - \L k (t )\ 2 L k (t ) 

+ L k {t )t k {t ) (L k {t )L k {t ) + L k {t )t k {t ) 
-- . 



(5.87) 
(5.88) 

(5.89) 
(5.90) 



5.5.1 Almost equilibrium states are Hadamard states 

Now we are prepared to prove that our construction indeed yields Hadamard states. 

Proposition 5.4 The two-point distribution of an almost equilibrium state upj satisfies 
the Hadamard condition. 

Proof. The strategy for the proof is the following. We show that that the difference of the 
two-point distributions, i.e., (u^ i — uUj)(x, x), is a smooth function and thus the singular 
parts of up j and u^j coincide. As indicated in section 12.2.41 an implicit infrared cutoff for 
the relevant integrals can be assumed, i.e., we have only to care for the large k behaviour. 
(For the closed case e = +1 no such assumption is needed.) Since the states of low energy 
u^j are known to have a Hadamard singularity structure, this proves the same for the 
almost equilibrium states upj. 

We have shown that bi = \/[S] and b 2 = 0. Hence, 



t ; (2) t ; (2) 
Upj ~ Uooj 



(x, x') 



dk F k (x)y k (x 



„ L k (x°)L k (x'°) + L k (x°)L k (x 

P /3ci(L k ,f) _ I 



I0\ 



(5.91) 



In |LR90j the following estimates for the growth of Y^(x) and T k (t) and their derivatives 
for large k are given (a is a multi- index, j e N): 

0(A;l Q l +2 ) , 



|Dx^k(x) 



and 



DlL k (t) = 0{k>-*) . 
Using this, we can then estimate the growth of the following expression 
sup |D^ [F k (x)F k (x') (L k (t)L k (f) +L k (t)L k (t'))] \ 
= sup |DWY k (x) D^' l F k (x / ) (D{L k (t) B{,L k (t') + DfZ k (t) D(,L k (t'] 

< |DL ff| ^k(x) 



D^V k (x') (\D 3 t L k (t)\ B{,L k (t') + \D{L k (t)\ D{,L k (t' 



(5.92) 

(5.93) 

(5.94) 
(5.95) 

(5.96) 
(5.97) 
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where \v\ = + \a'\ and a = \a\ +j, a' = \a'\ We also need an estimate for the growth 
of Ci(k) := Ci(T k , /). In [OlbOTb] it has been shown that the k-dependence of Ci(T k , /) for 
large k is such that there exist constants a-i > with 

ai(l + k) < ci(k) < a 2 (l + k) . (5.98) 

Thus, the factor -j^jz^j) ^ , which vanishes faster than any polynomial, makes the integral 

- uS,t)M (5.99) 

converge absolutely. This proves that the two-point difference ( 15 . 9 1 [) is a smooth function 
which in turn proves the Hadamard property of lo\ Jx, x'). □ 
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6 Summary and Outlook 



The present thesis is the successful accomplishment of a natural task suggested by recent 
results in quantum field theory in curved spacetimes. In |Olb07a] ground states, so-called 
states of low energy, for the worldline averaged, renormalized stress-energy tensor of the 
Klein-Gordon field in Robertson- Walker spacetimes were constructed. They are obtained 
by minimizing the averaged energy density that an isotropic observer measures in a pure, 
homogeneous, isotropic, quasi-free state. The states of low energy depend on the sampling 
function used in the averaging procedure. 

In this thesis we constructed, in the same setting, a family of states that we consider as 
almost equilibrium states to the inverse temperature (3. The almost equilibrium states are 
obtained by a suitable application of the principle of minimal free energy - a cornerstone 
of statistical mechanics. More precisely, we defined a free energy functional on the homo- 
geneous, isotropic, quasi-free states of the Klein-Gordon field, where the 'inner energy' was 
constructed from the worldline averaged stress-energy tensor since that quantity is known to 
be lower bounded on the class of Hadamard states by a quantum energy inequality |FewOO] . 
We showed how this energy can be written in terms of the same parameters that determine 
the two-point distribution of a homogeneous, isotropic, quasi-free state. Viewing each mode 
of the quantum field as a quantum system with one degree of freedom, we assigned to the 
modes the usual von Neumann entropy. An essential point of the whole construction was 
the determination of the entropy of such a system in terms of the two-point distribution. 
This was accomplished by calculating the generator of a KMS state on the Weyl algebra 
of a system with a single degree of freedom. By subsequent minimization of the free en- 
ergy, we obtained an explicit expression for the two-point distribution of a family of states, 
determined, again, by the sampling function of the averaging procedure and the inverse tem- 
perature j3. Finally, we showed that the two-point distributions of the almost equilibrium 
states satisfy the Hadamard condition. The latter step is a vital part of the construction 
as otherwise the utilized quantum inequality would not guarantee the existence of a finite 
lower bound for the averaged stress-energy. 

The formulation of quantum field theory in curved spacetimes disentangles the construc- 
tion of the quantum algebra from the construction of the quantum states. It is well known 
how to construct the algebra of a linear quantum field in a globally hyperbolic spacetime. 
The construction of physically meaningful states, however, is less straightforward. With 
decreasing symmetry properties of the underlying spacetime it becomes increasingly com- 
plicated to single out such states or rather classes of them. Regarding timelike symmetries, 
the step from stationary to non-stationary spacetimes is the biggest one. It brings along 
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unitarily inequivalent representations of the algebra of observables and makes the concept 
of a conserved Hamiltonian that generates global time evolution useless. In view of this 
situation, it is important to watch for other notions of energy that may be more meaningful. 
We advocate that such a notion is provided by the averaged energy densities used in the 
quantum energy inequalities. 

The almost equilibrium states defined here are, to date, the sole example of explicitly 
constructed (global) Hadamard states with thermal properties on a non-stationary space- 
time. For quantum fields in a non-stationary spacetime it is not reasonable to expect the 
existence of true equilibrium states with definite temperature 1/(3 since the influence of 
tidal forces will always destroy such property. We claim that the almost equilibrium states 
are reasonably defined approximations to equilibrium, and that they provide an interesting 
starting point for future investigations. Due to the setting, namely, quantum fields in cos- 
mological spacetimes, they might provide a class of states useful in the inflationary scenario 
and the analysis of the cosmic microwave background. 

Like the states of low energy, the two-point distribution of an almost equilibrium state 
involves a sampling function f(t) e @(R). The function f(t) describes the measurement of 
energy by an isotropic observer. This function introduces some freedom in the construction 
which may be used to design states with desired properties. This might be useful with 
respect to models of the early universe. 

Admittedly, our construction of almost equilibrium states makes heavy use of the mode 
decomposition, which is available in Robert son- Walker spacetimes because of the underlying 
symmetries. In general, there is no satisfactory mode decomposition on an arbitrary non- 
stationary spacetime, and thus there is no direct route for a generalization of our method. 

Quantum energy inequalities, which play a major role in our construction, are among 
a family of interconnected criteria for dynamical stability of quantum systems [FV03J. In 
[SVOH], a link between local thermal equilibrium states of a linear scalar field on a curved 
spacetime and quantum energy inequalities was alluded. It was proved that the existence of 
a linear scalar quantum field fulfilling some local thermal condition, in the sense of [BOR02J, 
implies a quantum energy inequality for these states. However, no states are known that 
fulfill the hypothesis of the paper. Due to our construction, which, starts with a quantum 
inequality rather than to derive one, it might turn out that the almost equilibrium states 
provide examples for local thermal equilibrium states (see below). 

Among the mentioned stability criteria for quantum systems there is passivity, which has 
not been generalized to non-stationary spacetimes so far. It is known that KMS states 
and mixtures of KMS states are passive. It might be worth considering, whether almost 
equilibrium states (and mixtures of almost equilibrium states) are 'almost passive' in a 
sense to be specified. Such a notion of passivity should be able to measure the amount of 
mechanical work that can be performed, owing to the present tidal forces, by a system in a 
cyclic process in a non-stationary spacetime. 
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Local thermal equilibrium states 

The compact support of the sampling function f(t) introduces intrinsically a kind of local- 
ization in time and, by taking the associated closed double cone, a localization in spacetime 
into the almost equilibrium states. An open question for future investigations might be how 
almost equilibrium states are connected to other notions of local thermal equilibrium. The 
most interesting states would be the local thermal equilibrium (LTE) states in the sense 
of Buchholz, Ojima, and Roos in [BOR02J. This approach gives a precise meaning to the 
saying that by a local measurement one cannot distinguish between local and global equi- 
librium states. This abstract idea is implemented in field theory by the following procedure. 
First, one chooses a set of global thermal equilibrium states, e.g., mixtures of KMS states 
in different inertial systems. Then, for a point x, one chooses a certain set S x of reference 
observables. A state uj is called ^-thermal if its expectations values for all observables 
in S x coincide with the expectation values for some global thermal equilibrium state uj x . 
Thus, as long as one considers observables from S x it is impossible to distinguish between 
the local state uj and the global thermal equilibrium state uj x . Of course, nothing is said or 
assumed about the expectation values of other observables. Thus, in general, uj does not 
coincide with the global state uj x . The spaces S x are chosen as the linear spaces generated 
by the balanced derivatives of the Wick-square : <fi 2 : (x). Balanced derivatives are defined 
in |BOR02j as 

&V : 2 : (x) = lim [<j>{x + ()<J)(x -()- cu o (0(x + Q^x - 0)1] , (6.1) 

where S M = d^n . . . d^ n with the multi-index /x := . . . , /x n ) e N n and uj is the vacuum 
state. Here, one takes the limit along spacelike directions Q G M in Minkowski spacetime, so 
that <p(x + C)4>i x ~~ C) is we U defined as a quadratic form [BosOOj . Now, with this approach 
one calculates for the Wick-square, which corresponds to the balanced derivative of zero-th 
order, of the massless field in Minkowski space 

^ e (: <f> : (*)) = " ( 6 - 2 ) 

independently of x, where e is an orthonormal tetrad defining the Lorentz frame of the 
KMS state (see the relativistic KMS condition [BB94J). This is interpreted as an indication 
that the the Wick-square serves as a scalar thermometer. 

It is natural to ask for the Wick-square of an almost equilibrium state and if there is 
a similar interpretation as thermometer. Hereby, one has to resolve an ambiguity in the 
definition the Wick-square since one has to choose a ground state. By construction, it is 
natural to choose the states of low energy for this purpose. But, by the method of locally 
covariant quantum fields [BFV03] , one could also choose the Minkowski vacuum. The latter 
procedure should be viable by the generalization of LTE states for curved spacetimes that 
is proposed [S V08j . 
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To the very end, let us make some speculative remarks regarding the former procedure. 
We use use the states of low energy for the parametrization and as the ground state. For 
simplicity, we consider flat Robertson- Walker spacetimes. By some formal manipulations, 
one easily, obtains the expression 

Urn ( u » - JZ) (« + C, x - = pip / <* ZtZ, ■ W 

Note that this expression bears a strong resemblance to its Minkowski spacetime counter- 
part. In that case by equation ( 13. 38ft we have 2|Lk(x°)| 2 = ^ and Ci(L k , f) = u; k provided 
that J f(t) 2 dt = 1. We learn from this that C\(T^ /) possibly plays the role of a general- 
ized frequency. In order to proceed into the direction of LTE states on needs estimates on 
Ci(L k , f). Such estimates may, e.g, lead to the interpretation that the almost equilibrium 
states are a mixture of LTE states with a range of temperatures a(f) < (3 < b(f). 
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